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Abstract  of  Dissertation  Presented  to  the  Graduate  School  of 
the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 

PHASE  SEQUENCE  OPTIMIZATION  IN  A NETWORK 
OF  MULTIPHASE  SIGNALS 

By 

Veretta  Johnson 
August,  1992 

Chairman:  Kenneth  G.  Courage 

Major  Department:  Civil  Engineering 

The  purpose  of  this  study  is  to  develop  a methodology 
which  will  determine  the  optimal  phase  sequences  at  each 
intersection  of  a network  of  multiphase  signals.  The  cycle 
length  and  green  time  allocation  are  assumed  to  be  already 
determined  and  are  not  optimized.  This  is  accomplished  in 
three  stages.  In  the  first  stage  an  algorithm  which 
determines  the  maximum  sum  of  inbound  and  outbound  green 
bands  on  a two  signal  arterial  highwway  is  developed.  This 
algorithm  gives  maximum  bandwidths  for  each  of  sixteen 
possible  phase  sequence  combinations  at  the  signals  of  the 
arterial  highway.  In  the  second  stage  the  algorithm  is 
extended  to  optimize  any  arterial  highway.  Optimal  phase 
sequence  combinations  and  corresponding  offsets  are  outputs 
of  the  algorithm.  Finally,  the  selection  of  the  best  phase 
sequence  combination  for  a closed  network  is  found  by  using 
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offsets,  splits,  and  cycle  length  to  formulate  a function  of 
several  variables  whose  maximum  value  is  the  maximum  sum  of 
green  bands  on  the  network.  A case  study  is  used  to 
demonstrate  the  selection  technique. 

The  results  of  this  study  may  be  used  to  supplement 
existing  traffic  signal  timing  plan  design  techniques  which 
do  not  include  phase  sequence  selection  capabilities. 
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CHAPTER  ONE 


INTRODUCTION 

Traffic  control  has  been  an  important  area  of  focus  for 
safety  on  roadways  since  the  start  of  the  automobile  era. 

The  earliest  control  devices  were  hand  signals  and  these 
signals  were  later  followed  by  traffic  lights.  While  the 
date  of  the  installation  of  the  first  electric  signal  in  the 
United  States  is  disputed  (Cleveland  in  19141,  Detroit, 
Michigan2,  Texas  in  19213)  , the  fact  that  signals  are  one  of 
the  most  fundamental  and  important  means  of  traffic  control 
is  not. 

The  evolution  of  traffic  signals  developed  in  two  major 
ways:  which  were  toward  coordination  and  in  the  direction  of 
traffic  responsiveness.  These  evolutionary  developments 
were  designed  to  accomplish  an  objective  other  than  the 
primary  one  of  safety.  Decreasing  the  inconvenience  of 
unnecessary  stops  and  delay  caused  by  traffic  signals  and 
facilitating  the  full  utilization  of  available  capacity  are 
two  of  their  objectives. 

Several  control  schemes  and  strategies  have  been 
developed  to  accomplish  these  objectives.  These  schemes 
include  traffic  signals,  traffic  management  measures,  one 
way  streets,  prohibition  and  route  control  devices.  The 
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control  strategies  include  the  principle  of  synchronization 
of  a group  of  traffic  signals.  Synchronization  is  achieved 
by  communication  between  individual  intersection 
controllers. 

The  maximum  through-band  synchronization  scheme  is  the 
oldest  and  most  widely  used  synchronization  scheme.  It 
consists  of  offsetting  the  beginning  of  the  green  phases  of 
successive  controllers  with  respect  to  a system  reference 
time  to  allow  as  many  vehicles  as  possible  to  drive  through 
these  intersections  at  some  design  speed  without  stopping. 
Bandwidth,  through-band  and  green-band  are  terms  used  to 
describe  the  number  of  seconds  in  each  cycle  in  which  a 
vehicle  traveling  the  arterial  street  is  able  to 
theoretically  start  at  the  first  signal  and  travel  through 
to  the  last  signal  at  the  design  speed  without  encountering 
a red  indication. 

The  time  space  diagram  is  a tool  used  for  the  analysis 
of  the  flow  of  vehicles  and  measurement  of  bandwidth  on  an 
arterial  street.  The  signals  of  the  arterial  highway  are 
represented  by  vertical  lines,  rectangles  superimposed  over 
the  vertical  lines  are  used  to  represent  red  (darkened 
rectangle)  and  green  (clear  rectangle)  periods,  and  lines  of 
constant  slope  are  used  to  represent  the  speed  of  the 
vehicles  alone  the  arterial.  One  axis  specifies  time  and 
the  other  axis  specifies  location  (space).  Figure  1.1 
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illustrates  a typical  time  space  diagram.  Each  direction  of 
an  arterial  highway  has  its  own  bandwidth. 
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INDICATES  OUTBOUND  TRAFFIC  FACES  RED  INDICATION 
INDICATES  INBOUND  TRAFFIC  FACES  RED  INDICATION 


FIGURE  1.1:  TYPICAL  TIME-SPACE  DIAGRAM 

Background  Discussion 

There  are  a number  of  bandwidth  maximization  methods. 
These  include  graphical,  algebraic  and  trial  and  error 
methods.  These  methods  were  limited  to  special  cases  such 
as  constant  block  lengths,  equal  green  bands  at  each 
intersection  and  constant  progression  speeds  if  the  system 
was  two-way.  Analytical  methods  have  been  developed  and 
programmed  which  have  more  flexibility.  Synchronization  is 
optimized  generally  by  maximizing  a weighted  sum  of  through- 
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bands  and/or  by  minimizing  a weighted  sum  of  stops  and 
delay.  The  timing  components  cycle  length,  splits,  offsets 
and  phase  sequence  are  adjusted  to  attain  the  optimum 
objective  function  value  in  both  synchronization  methods. 

Network  control  is  a form  of  signal  control  for  a group 
of  adjacent  signalized  intersections  where  coordination  is 
considered.  The  basic  terms  used  in  discussing  signalized 
intersection  control  are  1)  phase — the  portion  of  a signal 
cycle  allocated  to  any  single  combination  of  one  or  more 
traffic  movements  simultaneously  receiving  the  right  of  way 
during  one  or  more  intervals,  2)  sequence — a series  of 
successive  phase  assignments  during  the  cycle,  3)  interval- 
-a  discrete  portion  of  the  signal  cycle  during  which  the 
signal  indications  remain  unchanged,  4)  cycle  length — the 
time  required  for  one  complete  sequence  of  signal  intervals, 
5)  split — the  percentage  of  a cycle  length  allocated  to  each 
of  the  various  phases  in  a signal  cycle,  and  6)  offset — the 
time  difference  between  the  start  of  the  green  indication  at 
one  intersection  as  related  to  the  system  time  reference 
point.  Figure  1.2  illustrates  these  components  of  signal 
phasing. 

Although  several  signal  phasing  components  were  defined 
previously,  this  research  concentrates  on  optimizing  phase 
sequences.  Phasing  is  a technique  to  reduce  conflicts 
between  traffic  movements  at  signalized  intersections. 
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FIGURE  1.2:  SIGNAL  PHASING  COMPONENTS 

The  phasing  issue  primarily  relates  to  a left  turn  phasing 
problem.  Left  turn  treatment  alternatives  include 
permitted,  protected  and  permitted/protected.  Permitted 
left  turn  treatment  indicates  that  the  left  turning  traffic 
is  permitted  to  make  the  left  turn  whenever  the  driver 
judges  that  the  gap  in  opposing  traffic  is  sufficient. 
Protected  treatment  indicates  that  left  turn  traffic  must 
wait  for  a green  arrow  to  turn.  Permitted/protected 
treatment  on  an  approach  indicates  that  both  types  of  left 
turn  treatments  are  used.  The  problem  of  determining  the 
type  of  left  turn  treatment  needed  has  been  explored  and 
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warrants  have  been  developed  to  aid  the  traffic  engineer  in 
making  this  decision4. 

In  analyzing  the  problem  of  bandwidth  maximization  only 
the  two  opposing  approaches  are  considered  at  a given  time 
for  each  intersection.  This  is  possible  because  the 
movements  from  the  side  approaches  can  not  move 
simultaneously  with  the  movements  on  the  approaches  under 
consideration.  For  notational  purposes  one  of  the 
directions  on  the  arterial  highway  under  consideration  will 
be  called  the  inbound  direction  and  the  opposing  direction 
will  be  called  the  outbound  direction.  The  following  phase 
sequences  use  phase  names  to  represent  which  movements  have 
the  right  of  way  during  the  phase.  The  "0"  designation 
means  that  both  the  outbound  movements,  through  and  left 
turning,  have  the  green  indication.  "I"  refers  to  the  phase 
in  which  both  inbound  movements  have  the  green,  "T" 
represents  the  phase  in  which  both  through  movements, 
inbound  through  and  outbound  through,  have  the  right  of  way 
and  "L"  refers  to  the  phase  in  which  both  left  movements 
have  the  green  indication.  Figure  1.3  illustrates  the 
movements  allowed  with  each  phase  name  described  above. 

The  following  phase  sequences  are  possible:  The  first 

half  of  the  name  of  each  sequence  reflects  whether  the 
outbound  left  turn  precedes  or  follows  the  through  movement 
and  the  second  refers  to  the  inbound  left  turn.  The  four 
phase  sequences  to  be  analyzed  are  as  follows: 
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FIGURE  1.3:  MOVEMENTS  ALLOWED  WITH  PHASES 

The  Lead-lag  Sequence 

In  this  sequence  the  outbound  left  turn  leads  and  the 
inbound  left  turn  lags.  The  phase  sequence  is  OTI  which 
means  outbound  left  and  through  have  the  green  first, 
followed  by  a second  phase  where  through  movements  have  the 
green  indication  and  in  the  final  phase  the  inbound  left  and 
through  movements  get  the  right  of  way.  Figure  1.4 
illustrates  this  phase  sequence. 

The  Lead-lead  Sequence 

Both  the  inbound  and  outbound  left  turns  lead  with  this 
phase  sequence.  This  phase  sequence  allows  one  of  the 
directions  to  be  assigned  more  time  than  the  other. 
Consequently  the  phase  sequence  is  either  LIT  or  LOT.  This 
phase  sequence  is  diagrammed  in  figure  1.5. 
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FIGURE  1.4:  THE  LEAD-LAG  SEQUENCE 
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FIGURE  1.5:  THE  LEAD-LEAD  SEQUENCE 
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The  Lag-lead  Sequence 

The  outbound  left  turn  lags  and  the  inbound  left  turn 
leads  with  the  lag-lead  sequence.  In  this  case  the  phase 
sequence  is  ITO.  Figure  1.6  illustrates  this  phase 
sequence. 

The  Lag-lag  Sequence 

Both  left  turns  lag  with  lag-lag  phasing.  The  phase 
sequence  is  TOL  or  TIL  depending  on  which  direction  requires 
more  time.  This  phase  sequence  is  illustrated  in  figure 
1.7. 

Statement  of  the  Problem 

Various  algorithms  have  been  developed  to  optimize  the 
timing  of  signalized  traffic  systems.  Cycle  length,  splits 
and  offset  optimization  algorithms  for  networks  have  been 
developed  and  have  been  programmed.  Algorithms  which 
completely  optimize  the  timing  components  of  arterial 
highways  have  also  been  developed.  The  existing  algorithms 
for  the  determination  of  phase  sequences  in  a network  of 
multiphase  signals,  however,  are  limited. 

Morgan  and  Little5  developed  a method  of  optimizing 
offsets  on  an  arterial  highway  using  the  cycle  length, 
splits  and  travel  times  on  the  arterial  highway.  Bleyl6 
developed  a model  for  selecting  cycle  length (within  ranges) 
and  offsets  which  maximize  the  bandwidth  on  an  arterial 
highway  using  travel  times,  volume  and  minimum  green  times. 
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Yardeni7,8  developed  a model  which  determines  cycle  length, 
splits  and  offsets  on  arterial  highways.  Brooks9  developed 
a method  which  selects  offsets,  splits  and  cycle  length  at 
signals  using  the  intersection  with  the  shortest  main  street 
green  (called  the  critical  intersection)  and  a range  of 
cycle  lengths.  Messer  et  al.10  developed  a computer  program 
for  Brook's  algorithm.  Hillier11  developed  an  algorithm  and 
Wagner  and  Gerlough12  developed  a computer  program  to  select 
offsets  on  an  arterial  highway.  Little  et  al13  developed  a 
model  which  optimizes  phase  sequences,  cycle  length  and 
splits  on  arterial  highways  and  triangular  networks. 

The  algorithms  for  analyzing  networks  are  much  more 
limited  and  only  one  optimizes  phase  sequence  combinations 
at  the  signals  of  the  network.  Robertson14  developed  an 
algorithm  which  determines  optimal  splits,  and  offsets  in  a 
multiphase  network  of  signals.  The  algorithm  has  been 
programmed  and  is  widely  used  in  optimizing  networks.  While 
an  optimal  cycle  length  can  be  determine  by  multiple  runs  of 
the  program,  phase  sequences  must  be  input  to  the  program 
and  as  such  are  not  selected  by  the  algorithm.  Messer  et 
al.15  developed  an  algorithm  which  optimizes  left  turn  phase 
sequences,  cycle  length,  and  offsets  in  signalized  networks. 
The  algorithm  requires  extensive  data  preparation  and 
substantial  computer  resources.16  Cohen17  and  Rogness18  have 
shown  that  optimization  of  left  turn  phasing  can 
significantly  improve  performance  on  signalized  networks. 
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Thus  this  paper  concentrates  on  the  need  to  develop  a method 
to  determine  the  best  phase  sequence  at  each  intersection  of 
a multiphase  network. 

In  this  research  the  goal  is  to  develop  an  algorithm 
which  determines  optimal  phase  sequences  for  the  signals  of 
a closed  network  which  can  be  implemented  using  limited 
computer  resources  and  which  when  combined  with  existing 
algorithms  provides  a complete  optimization  of  a network  of 
multiphase  signals. 

The  problem,  therefore,  is  to  assign  appropriate  phase 
sequences  to  the  signals  of  a multiphase  network  so  that  the 
overall  inconvenience  of  unnecessary  stops  and  delay  is 
minimized  and  progression  on  the  arteries  of  the  network  is 
maximized.  The  objective  of  this  research  is  to  develop  an 
algorithm  which  selects  optimal  phase  sequences  at  each 
intersection  of  a network  of  multiphase  signals  using  the 
maximum  sum  of  green  band  as  the  performance  criterion  and 
which  requires  minimal  data  input  and  computer  resources. 

Scope  of  the  Study 

This  dissertation  focuses  on  developing  an  algorithm 
which  determines  the  optimal  phase  sequence  combinations  in 
a network  of  multiphase  signals.  The  study  does  not  attempt 
to  present  guidelines  for  selecting  the  appropriate  type  of 
left  turn  signal  treatment  or  to  optimize  cycle  length, 
green  allocation  or  offsets. 
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Figure  1.8  depicts  the  analytical  steps  carried  out  in 
the  design  of  a coordinated  network  of  multiphase  signals. 
Even  though  one  can  argue  that  the  determination  of  green 
apportionment (splits)  could  be  coincidental  with  phase 
sequence  selection,  it  is  assumed  in  this  research  that 
splits  are  known. 

Gartner19-20  has  found  that  since  both  throughput  capacity 
and  delay  time  increase  with  cycle  length,  the  common  cycle 
time  should  be  the  minimum  required  to  satisfy  the 
intersection  with  the  largest  traffic  volume.  Other  methods 
of  determining  network  cycle  time  have  been  described  by 
Webster21  and  in  the  Highway  Capacity  Manual4.  For  purposes 
of  this  research  network  cycle  time  is  assumed  to  be 
determined. 

The  linking  of  signals  is  obtained  by  appropriate 
selection  of  offsets  through  the  network.  The  network 
closure  constraint  requires  that  the  sum  of  the  offsets 
around  each  loop  of  the  network  be  an  integer  multiple  of 
the  cycle  length.  Network  optimization  programs 
TRANSYT-7F22,  MAXBAND-8  615  and  SIGOP  II23-24  all  provide  good 
solutions  to  the  offset  optimization  problem.  Offset 
optimization  is  not  within  the  scope  of  this  research 


effort. 
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FIGURE  1.8:  STEPS  IN  THE  DESIGN  OF  A COORDINATED  NETWORK 

OF  MULTIPHASE  SIGNALS 

To  accomplish  the  goals  of  this  research,  the  following 
analytical  steps  were  taken. 

1)  Develop  a model  which,  given  the  green  requirement 
per  movement,  determines  the  green  time  per  phase 
for  each  phase  sequence. 

2)  Develop  a method  of  determining  the  beginning  and 
end  of  green  at  an  intersection  as  a function  of 
the  offset  at  the  intersection  for  each  phase 
sequence. 

3)  Develop  a closed  form  expression  for  the  outbound 
and  inbound  bandwidth  between  pairs  of  consecutive 
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intersections  as  a function  of  the  offset  at  the 
intersections . 

4)  Develop  a method  of  determining  the  offset  values 
which  yield  maximum  value  of  the  sum  of  the 
bandwidth  functions  developed  under  objective  3. 

5)  Develop  an  algorithm  to  determine  the  phase 
sequence  pattern  which  yields  the  maximum 
attainable  green-band  on  an  artery  using  the  values 
obtained  in  objective  4. 

6)  Develop  a method  of  putting  solutions  to  the 
arteries  of  a network  together  so  as  to  optimize 
the  sequences  for  the  entire  network. 

These  objectives  require  explicit  consideration  of 
several  complex  coordination  problems. 

Organization 

In  chapter  two  a step  by  step  procedure  for  determining 
the  maximum  sum  of  inbound  and  outbound  bands  on  a two 
intersection  arterial  highway  is  developed.  In  chapter 
three  a discussion  of  a method  of  successively  treating  the 
constraints  imposed  by  each  signal  of  the  arterial  highway 
which  allows  the  extension  of  the  two  signal  arterial 
highway  solution  to  any  arterial  highway  is  given.  Chapter 
four  provides  a means  of  extending  the  arterial  solution  to 
a general  closed  network.  Conclusions  and  recommendations 
are  given  in  chapter  five. 


CHAPTER  TWO 


DEVELOPMENT  OF  AN  ALGORITHM  FOR  THE 
TWO  SIGNAL  ARTERIAL  HIGHWAY 

In  this  chapter  an  algorithm  to  determine  the  optimal 
solution  to  systems  consisting  of  two  signals  is  developed. 
The  first  step  is  to  determine  the  number  of  seconds  of 
green  to  assign  to  each  phase.  The  green  per  phase 
assignments  are  then  used  to  determine  the  beginning  and  end 
of  the  green  time  at  each  intersection.  Using  these  values 
the  bandwidth  between  two  signals  is  expressed  as  a function 
of  the  offset  at  the  second  signal  of  the  system.  Finally, 
a closed  form  solution  to  the  problem  of  finding  the  maximum 
sum  of  the  outbound  and  inbound  bandwidths  is  developed. 

Determination  of  Green  Time  Per  Phase 
at  a Single  Intersection 

In  general,  coordinated  systems  operate  best  when  the  main 
street  flow  is  predominantly  through  traffic.  Therefore  to 
achieve  the  goal  of  maximizing  the  sum  of  the  inbound  and 
outbound  bandwidths  the  objective  in  assigning  green  time  to 
the  phases  is  to  assign  as  much  time  to  the  through  movements 
as  possible.  The  number  of  seconds  of  green  to  be  allocated 
to  each  phase  can  be  determined  using  the  green  time 
requirement  per  movement.  The  movements  are  "OBL", 
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representing  the  outbound  left  turn  green  time  requirement, 
"OBT"  representing  the  outbound  through  movement  green  time, 
"IBL"  representing  inbound  left  turn  green  and  "IBT" 
representing  the  inbound  through  green  time  requirement. 
These  representations  allow  the  following  assignments  to  be 
made  for  each  of  the  four  phase  sequences  previously 
considered. 

1)  For  lead-lag  phasing,  the  phase  sequence  is 
identified  as  OTI  which  requires  an  assigned  time  to  the 
outbound  movement  0,  the  through  movement  T,  and  the  inbound 
movement  I.  With  this  phase  sequence  the  outbound  left 
traffic  is  served  only  during  the  0 phase.  Therefore  the 
number  of  seconds  assigned  to  the  0 phase  must  be  at  least 
OBL.  In  order  to  achieve  the  goal  of  assigning  as  much  time 
to  the  through  movement  as  possible,  the  decision  made  here 
is  to  limit  the  duration  of  the  0 phase  to  be  exactly  OBL. 
Similarly  since  the  inbound  left  turn  traffic  is  only  served 
during  the  I phase,  the  decision  is  to  assign  a number  of 
seconds  at  least  equal  to  IBL  to  the  I phase  and  to  limit 
the  duration  so  that  it  is  exactly  equal  to  IBL.  Hence,  the 
assignment  of  green  time  to  the  phases  for  lead-  lag  phasing 
are  as  follows: 

0 = OBL  = number  of  seconds  of  green  required  for  the 

outbound  left  turning  traffic 

1 = IBL  = number  of  seconds  of  green  required  for  the 

inbound  left  turning  traffic 

T = Cm  - OBL  - IBL 
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If  T < 0 then  let  0=0,  1=0,  and  T = Cm 
where  Cm  is  the  number  of  total  seconds  allocated  to  the 
outbound  and  inbound  movements  on  the  artery. (Note  C - Cm  is 
the  total  number  of  seconds  allocated  to  the  cross  street 
movements  at  the  intersection.)  The  "If"  statement  causes 
the  left  turn  protection  to  be  eliminated  in  the  event  that 
the  total  number  of  seconds  allocated  to  the  artery  is 
insufficient. 

2)  For  lag-lead  phasing  the  assignments  of  green  time 
per  phase  are  made  the  same  as  for  lead-lag  phasing.  3) 

For  lead-lead  and  4)  lag  -lag  phasing  the  phases  used  are 
L,0, I,  and  T.  The  assignment  of  green  time  to  each  phase  is 
made  with  the  goal  of  maximizing  the  number  of  seconds 
assigned  to  the  through  movement.  The  number  of  seconds 
assigned  to  each  phase  is 
L = min  [OBL,  IBL] 

I = max  [0,  IBL  - L] 

0 = max  [0,  OBL  - L] 

T = Cm-L-I-0 

If  T < 0 then  set  1=0,  0=0,  L=0  and  T = Cm. 

Determining  the  Beginning  and  End  of  Green 
The  key  requirements  in  calculating  the  bandwidth  on  an 
arterial  highway  are  the  determination  of  the  beginning  and 
end  of  the  green  period  at  each  intersection.  The  beginning 
and  end  of  green  at  a signal  depend  on  the  offset,  phase 
sequence  and  green  time  per  phase. 
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Let  a0  represent  the  beginning  of  the  outbound  green, 
b0  the  end  of  the  outbound  green, 
a;  the  beginning  of  the  inbound  green,  and 
bj  the  end  of  the  inbound  green. 

These  representations  determine  the  beginning  and  end  of  green 
in  each  phase  sequence  in  the  following  ways: 

1)  Lead-Lag  Phasing 

Figure  2.1  illustrates  the  phases  of  the  lead-lag  phase 
sequence.  The  inbound  and  outbound  directions  are  labeled  on 
the  right  of  the  figure.  The  outbound  green  begins  at  the 
offset  and  continues  until  the  time  that  the  inbound  left  turn 
traffic  is  given  the  green  indication.  The  I phase  begins  at 
this  time  and  is  labeled  b0  on  the  figure.  The  inbound  green 
begins  at  the  time  labeled  a(.  The  inbound  time  represents  the 
time  that  the  outbound  left  turning  traffic  is  stopped  and  the 
inbound  through  traffic  is  then  given  the  green  indication. 
The  inbound  green  continues  through  the  I phase  and  ends  at 
the  point  labeled  bj  where  the  cross  street  traffic  is  given 
the  green  and  all  traffic  on  the  arterial  highway  is  stopped. 

For  the  lead-lag  sequence  a0,  b0,  ait  and  bj  are  computed  as 
follows : 

a0  = offset 

bQ  = offset  + O + T 

aj  = offset  + 0 

bj  = 


offset  + 0 + I + T 
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note:  indicates  that  the  through  traffic  faces  a red 

indication 

FIGURE  2.1:  BEGINNING  AND  END  OF  GREEN  FOR 

LEAD-LAG  PHASE  SEQUENCE 

The  values  of  a0/  b0,aj,  and  bj  are  determined  similarly  for 
the  remaining  three  phase  sequences  using  figures  2.2,  2.3  and 
2.4.  The  resulting  values  are  as  follows: 

2)  Lead-Lead  Phasing 

With  this  phase  sequence  either  0=0  or  1=0.  In  either 
situation  the  expressions  below  will  give  the  correct  values 
for  each  of  the  phases  a0,  b0,  bj,  and  bj  as  illustrated  in 
figure  2.2. 

a0=  offset  + L + I 

b0  = offset  + L + I + 0 + T 

a;  = offset  + L + O 

bj  = offset  + L + 0 + I + T 
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OR 


a; 


inbound 


■T- 


I 


outbound 


offset 


a, 


note:  indicates  that  the  through  traffic  faces  a red 

indication 


FIGURE  2.2:  BEGINNING  AND  END  OF  GREEN  FOR 

LEAD-LEAD  PHASE  SEQUENCE 


3 ) Lag-Lead  Phasing 
a0=  offset  + I 

bQ  = offset  + I + 0 + T 

a;  = offset 

bj  = offset  + I + T 

4 ) Lag-Lag  Phasing 
a0=  offset 

bQ  = offset  + T+  0 
a;  = offset 


b;  = offset  + T + I 
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*>i 


inbound 


outbound 


offset  a0  b0 

note:  |^H  indicates  that  the  through  traffic  faces  a red 

indication 


FIGURE  2.3:  BEGINNING  AND  END  OF  GREEN  FOR 

LAG-LEAD  PHASE  SEQUENCE 


inbound 


outbound 


OR 


b; 


inbound 


| outbound 


indicates  that  the  through  traffic  faces  a red 
indication 


FIGURE  2.4:  BEGINNING  AND  END  OF  GREEN  FOR 

LAG-LAG  PHASE  SEQUENCE 
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Determining  the  Bandwidth  Between  Two  Signals 

The  bandwidth  between  two  signals  can  be  computed  from 
the  beginning  and  end  of  the  green  at  each  of  the  two 
signals. 

Let  a(0l)  represent  the  beginning  of  the  outbound  green 
at  signal  1, 

b(0,i)  the  end  of  the  outbound  green  at  signal  1 
a(U)  the  beginning  of  the  inbound  green  at  signal  1 
b(jl)  the  end  of  the  inbound  green  at  signal  1 and 
N the  offset  at  the  second  signal. 

(a(o2),  b(o2)/  a(i2)/  and  b(i2)  are  similarly  defined  respectively 
for  signal  2) . Also  let  TTO  represent  the  outbound  travel 
time  between  signals  1 and  2 and  let  TTI  represent  the 
inbound  travel  time  between  the  two  signals.  Figure  2.5 
depicts  how  these  values  together  with  the  travel  time  can 
be  used  to  determine  the  bandwidth. 

Figure  2.5  shows  that  the  bandwidth  is  given  by 
EOB  - BOB,  where  EOB  = end  of  band  = min  [b(0l)/b(02)  -TT]  and  BOB 
= beginning  of  band  = max  [a(0l),  a(o2)  - TT]  . 

The  outbound  and  inbound  bandwidth  between  any  pair  of  signals 
can  be  computed  from  the  beginning  and  end  of  green  at  each  of 
the  signals  along  with  the  travel  time  between  the  two 
signals . 

The  next  step  is  to  examine  each  of  the  possible  sixteen 
combinations  of  phase  sequences  at  the  two  signals, 
examination  gives  the  following  outbound  and  inbound 


This 
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DISTANCE  (SPACE) 


GREEN  BAND  - MIN  { b,  + TT,  b2  } - MAX  { a,  + IT,  a*  } 


FIGURE  2.5  TIME  SPACE  DIAGRAM  SHOWING  THE  BANDWIDTH  BETWEEN 

TWO  SIGNALS 

calculations.  For  convenience  let  the  offset  at  the  first 
signal  be  zero.  Then  the  following  values  are  defined  for 
combinations  of  phase  sequences  at  the  two  intersections. 
Case  1)  Lead-Lag  at  Signal  1 / Lead-Lag  at  Signal  2 
AT  SIGNAL  1 AT  SIGNAL  2 


l(0,l) 

= 0 

a(o,2)  = 

N 

’(0,1) 

= o, 

+ 

T, 

b(o,2)  = 

N 

+ 

02 

+ 

T2 

(i.l) 

= 0, 

a(i,2)  = 

N 

+ 

°2 

'(i.l) 

= T, 

+ 

Ii  + °i 

II 

CS 

N 

+ 

T2 

+ 

i2  +o2 

The  outbound  beginning  of  the  green  band  denoted  here  as  B0 


B0  = max  [OB,,OB2  -TTO] 
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max  [0,  N-TTO] 


The  outbound  end  of  the  green  band  denoted  here  as  E0  is 


If  0,  + T,  > 02  + T2,  then  the  numbers  TTO  and  TT0+0,+T,-  02-T2 
divide  the  number  line  into  three  parts  as  shown  in  figure 
2.6.  The  values  of  B0  and  E0  are  determined  for  values  of 
the  offset,  N,  within  each  section  of  the  number  line  as 
indicated  in  figure  2.6. 

Using  the  expressions  for  B0  and  E0  one  may  then  obtain 
expressions  for  the  outbound  bandwidth,  OBW.  The 
expressions  are  shown  in  figure  2.7  for  each  value  of  N,  the 
offset  at  signal  2.  Therefore  the  outbound  bandwidth  can  be 
expressed  as  a function  of  N as 


If  0,+T,  < 02+T2,  then  the  number  line  is  again  divided 
into  three  sections  as  shown  in  figure  2.8  and  the  values  of 
B0  and  E0  are  as  shown  in  the  figure.  Using  these  values  the 
value  of  the  outbound  bandwidth  in  this  case  is  illustrated 
in  figure  2.9  and  can  be  written  in  function  form  as 


E0  = min  [ b(0>1) , b(o2)  - TTO] 

= min  [Oj+T,,  N+02+T2-TT0] 


OBW (N) 


N + 02  + T, 
02  + T2 


TTO  N < TTO 

TT0<  N < TT0+O!+T,-O2-T2 
N > TT0+01+T,-02-T2 . 


O.+T.+TTO-N 


OBW (N) 


N + 02  + T2  - 
0,  + T, 


0,  + T,  + TTO 


TTO 


N 


N > TTO 


N < TT0+0,+T,-02-T2 
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Now,  let  A = min  [ TTO  , TT0+0,+T,-02-T2  ] 

B = max  [ TTO  , TTO+Oj+Tj-Oj-Tj  ] 
and  C = 02  + T2  - TTO. 

Then  the  outbound  bandwidth  function  OBW  is  defined  as 
follows : 


OBW (N) 


N+C 

A+C 

A+B+C-N 


N < A 
A < N < B 
N > B. 


The  inbound  bandwidth  function  is  determined  using  a 
similar  analysis  of  the  values  of  N and  is  as  follows: 


If  Ii+T,  > I2+T2,  then 


IBW(N) 


' N+I2+T2-TTI+02-0,  N < 

• I2+T2  * TT+0,-02<  n < 

o,+i1+t,+tti-o2-n  n > 


tti+o,-o2 

TTI+O^Oj+I^T,-^-^ 

TTI+Oj-Oj+Ij+Tj-Ij-^ 


If  1,+Tj 
IBW(N)  = - 


< I2+T2,  then 

N+I2+T2-TTI+02-0!  N < TTI+0i-02+I,+T,-I2-T2 

I,+T,  TTI+O^O^I^T,-^-^  N < TTI+0,-02 

O^Ij+T^TTI-Oj-N  N > TTI+O^Oj 


Now  if  A,  B,  and  C are  defined  as 

A = min  [TTI+0,-02,  TTI+0,-02+I1+Ti-I2-T2] 

B = max  [TTI+0,-02,  TTI  +O,-02+I1+T1-I2-T2] 
C = 02  + I2  + T2  - 0,  - TTI 


again  the  bandwidth  function  can  be  expressed  as 


N < A 
A < N < B 
N > B. 


IBW(N) 


N+C 

A+C 

A+B+C-N 
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B0{ 

1_ 

0 N - 

— 1 

TTO 

N - 

TTO 

1 

E0{ 

0 

N+02+T2-TT0  1 N+02+T2-TT0  Oj+Tj 

TTO  TT0+(O,+T1-O2-T2) 

h 

a)  IF 

TTO  < TTO+  (01+T1-02-T2) 

OR 

B0 

h- 

0 

— 1 

0 

N 

- TTO 

I 

E0 

1 

N+02+T2-TT0  Oi+T, 

0 TT0+(01+T1-02-T2) 

0,+T, 

TTO 

\~ 

b)  IF 

TTO  > TT0+(0,+T1-02-T2) 

FIGURE  2.6:  EXPRESSIONS  FOR  THE 

OUTBOUND  BAND  LEAD-LAG 

BEGINNING  AND  END 
PHASE  SEQUENCE 

OF  THE 

OBW 

1 

N+02+T2-TT0  ( 02+T2  Oi+T(+TTO— n 

1 

0 

H 1 

TTO  TT0+(0i+T,-02-T2) 

a)  IF  TTO  < TT0+(0,+Ti-02-T2) 

OBW 

1 

N+02+T2-TT0  0,+T,  0,+Ti  +TTO— N 

1 1—  

TT0+(0,+T,-02-T2)  TTO 

b)  IF  TTO  > TT0+(0,+T,-02-T2) 


FIGURE  2.7:  VALUE  OF  THE  OUTBOUND  BANDWIDTH  FOR  LEAD-LAG 

PHASE  SEQUENCE 

Thus  both  the  outbound  and  inbound  bandwidth  functions  fit  a 
general  form  as  illustrated  above.  The  next  case  to  be 
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examined  is  one  where  the  phase  sequence  at  signal  1 is 
lead-lag  and  the  phase  sequence  at  signal  2 is  lead-lead. 
The  bandwidth  functions  are  calculated  as  shown  below  and 
fit  into  the  same  general  form  as  the  preceding  case. 

Case  2)  Lead-Laa/Lead-Lead  Phasing 
AT  SIGNAL  1 AT  SIGNAL  2 

a(o,2)  = N + Lj  + I2 
k(o,2)  = N + L>2  +I2+02+T2 
a(i,2)  = N + Lj  + 02 
b(j ,)  = Tt  + Ij  b(!2)  = N + L>2  +02+I2+T2 


a(o,D  = 0 


^(o.i)  - + T, 


a(i,D  = 


The  outbound  beginning  of  band,  B0  = max  [0B1,0B2  - TTO] 

= max  [0,  N+L2+I2-TTO]  . 

The  outbound  end  of  band,  E0  = min  [b(0l),  b(o2)-TTO] 

= min  [Oj+Tj,  N+L2+I2+02+T2-TT0]  . 
Thus  if  0,  + T,  > 02  + T2,  then  the  values  of  Bc  and  E0  are  as 
indicated  in  figure  2.8  and  the  value  of  the  outbound 
bandwidth  is  as  in  figure  2.9.  The  functional 
representation  of  the  bandwidths  are  as  follows: 


OBW(N)  = 


N+L2+I2+O2+T2-TT0  N < TTO-Lj-Ij 

02  + T2  TTO-L2-I2  < N < TT0-L2-I2+0,+T,-02-T2 

o,+t1+tto-l2-i2-  n n > tto-l2-i2+o1+t,-o2-t2 


and  if  0,+T,  < 02  + T2,  then 
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OBW(N) 


N+L2+I2+02+T2-TT0  N < TT0-L2-I2+01+T,-02-T2 

O,  + T,  ’ ‘ TT0-L2-I2+0i+T,-02-T2  < N < TTO-Lj-Ij 

0,+T,+TT0-L2-I2-  N N > TTO-L2-I2 


The  inbound  bandwidth  IBW  is  computed  as  follows: 


If  0,+1,+T,  > I2+T2,  then 


IBW (N) = • 


N+Lj+Oj+Ij+Tj-TTI-O, 

I2+T2  TTI+O, -I^-Oj  < 

0,+I1+T,-L2-02+TTI-N 


N < TTI  +0,-L2-02 
N < TTI+O,  -Lj-Oj+I  , +T,  - I2~T2 
N > TTI+01-L2-02+I1+T,-I2-T2. 


If  I,+T,  < I2+T2,  then 


IBW (N) 


N+L2+02+I2+Ti-TTI-01  N<  tti+o,-l2-o2+i,+t,-i2-t2 

I,+T,  ‘ ' 'TTI+0,+I,+T,-I2-T2  < N < TTI+O, 

0,+I,+T,+TTI-N  N > TTI+O,. 


B0  o n+o2+t2-tto  n+o2+t2-tto 

1 1 

E0  N+L2+I2+02+T2-TT0  N+L2+I2+02+T2-TT0  o,+t, 

TTO-Lj-Ij  TTO-Lj-Ij+COj+Tj-Oj-Tj) 


FIGURE  2.8:  EXPRESSIONS  FOR  THE  BEGINNING  AND  END  OF 

OUTBOUND  BANDWIDTH  FOR  LEAD-LEAD  PHASE  SEQUENCE 


OBW  N+L2+I2+02+T2-TT0  02+T2  0,+T,+TT0-L2-I2 


0 


TTO-I>2-I2  tto-l2-i2+o,+t,-o2-t2 


FIGURE  2.9:  VALUE  OF  THE  OUTBOUND  BANDWIDTH  FOR  LEAD-LEAD 

PHASE  SEQUENCE 
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General  Form  of  the  Bandwidth  Function 


The  outbound  and  inbound  bandwidth  functions  are  all  of 
the  general  form  discussed  below. 

Let  Fj(x)  be  a piecewise  linear  function  defined  by 


The  maximum  value  of  this  type  function  is  easily  seen 
to  be  attained  at  any  x in  the  interval  [A^BJ  and  is  equal 
to  Aj+C;.  Figure  2.12  illustrates  what  the  graph  of  a 
general  bandwidth  function  looks  like. 

The  solution  to  the  maximum  sum  of  two  such  functions 
depends  on  the  relationship  of  the  A/  s to  the  B/s.  There 
are  six  cases  to  consider  ( since  A;  < Bj)  . 

CASE  1.  A,  < B,  < A2  < B2 

CASE  2 . A,  < A2  < B,  < B2 

CASE  3.  At  < A2  < B2  < B, 

CASE  4.  A2  < B2  < A,  < B, 

CASE  5.  A2  < A,  < B2  < Bj 

CASE  6.  A2  < A,  < B,  < B2 

The  maximum  value  of  the  sum  of  F,  and  F2can  then  be 
determined  for  each  case.  Figure  2.13  illustrates  the 
solution  for  case  1.  The  sum  , F,+F2,  is  defined  in  this 
case  as  follows: 


F;  (X)  = 


X + C; 

a,  + q 

A;  + C;  + Bj  - X 


X < Aj 
A;  < X < Bj 
X > B;. 


2X  + C,  + C. 


1 T '-2 


x < A, 


X + A,  + C,  + C2 
A[  + B[  + Cj  + C2 
A,  + B,  + C!  + A2  + C2  -X 
Aj+Bj  +C,  + A2  + B2  + C2-2X 


A,  < X < B 

B,  < X < A2 
A2  < X < b2 


X >=B2 . 


(F,+F2)  (X) 
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The  maximum  value  of  the  sum  is  given  at  X in  the  interval 
[B,,A2]  with  the  value  of  the  sum  equal  to  Aj+Bj+Cj+Cj.  The 
maximum  for  each  case  described  above  is  shown  in  the  table 
below. 


Case 


Maximum  Sum  optimal  X 


1.  A, 

< 

B, 

< 

a2 

< 

b2 

C, 

+ 

c2 

+ 

a, 

+ 

B, 

B, 

< 

X 

< 

A 

2.  A, 

< 

a2 

< 

B, 

< 

b2 

C, 

+ 

c2 

+ 

A, 

+ 

a2 

a2 

< 

X 

< 

B 

3.  A, 

< 

a2 

< 

b2 

< 

B, 

C, 

+ 

c2 

+ 

A, 

+ 

a2 

a2 

< 

X 

< 

B 

4.  A-, 

< 

b2 

< 

A, 

< 

B, 

C, 

+ 

c2 

+ 

a2 

+ 

b2 

b2 

< 

X 

< 

A 

5.  A2 

< 

A, 

< 

b2 

< 

B, 

C, 

+ 

C2 

+ 

A, 

+ 

a2 

A, 

< 

X 

< 

B 

6 • A2 

< 

A, 

< 

B, 

< 

b2 

C, 

+ 

c2 

+ 

A, 

+ 

a2 

A, 

< 

X 

< 

B 

2 

1 

2 
1 
2 
1 


For  each  of  the  sixteen  phase  sequence  combinations  the 
maximum  sum  of  inbound  and  outbound  bandwidths  is  computed 
using  the  solution  for  the  general  function  form  and  the 
definitions  of  the  constants  A;,  B;,  and  C{  for  each  of  the 
phase  sequence  combinations.  The  result  is  that  the  maximum 
sum  of  the  inbound  and  outbound  bands  is  determined  for  each 
of  the  sixteen  phase  sequence  combinations. 

Thus  the  maximum  sum  of  inbound  and  outbound  bandwidth 
between  any  two  signals  can  be  determined  by  first  defining 
the  outbound  bandwidth  function  and  the  inbound  bandwidth 
function.  After  defining  the  two  functions  the  maximum  sum 
can  be  found  as  demonstrated  above. 
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FIGURE  2.10:  GENERAL  BANDWIDTH  FUNCTION 


FIGURE  2.11 


SUM  OF  TWO  BANDWIDTH  FUNCTIONS 
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With  phase  sequence  combination  1)  lag-lag/ lag-lag,  for 
example,  A,  = TTO, 

B,  = TT0+T,+0!-T2-02, 

C,  = T2+02-TT0 

A2  = TTI 

B2  = T,+Ij+TTI-T2-I2 

Cj  = t2+i2-tti. 

The  expressions  for  the  inbound  and  outbound  bandwidth 
for  each  of  the  remaining  possible  phase  sequence 
combinations  are  derived  similarly.  The  expressions  for 
these  functions  depend  on  the  phase  sequences  at  each  end  of 
the  artery.  Values  for  A;,  Bj,  and  Ct  for  each  of  the  sixteen 
possible  phase  sequences  are  shown  in  appendix  A. 


SUMMARY  OF  AN  ALGORITHM  OPTIMIZING  A TWO  SIGNAL  SYSTEM 

The  optimal  phase  sequence  combinations  for  each  signal 
of  a two-signal  artery  can  be  summarized  in  the  following 
steps. 

Step  1.  Do  the  following  at  each  signal 
Do  the  following  for  each  of  four  phase  sequences 
(lead-lag,  lead-lead,  lag-lead,  and  lag-lag) . 

a.  Assign  green  time  to  each  phase  using  formulas 
2.1,  2.2,  2.3,  and  2.4. 

b.  Determine  the  beginning  and  end  of  the  green  as 
functions  of  the  offset  at  signal  2 using  formulas 
2.5,  2.6,  2.7  and  2.8. 
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Step  2.  Do  the  following  for  each  combination  of  phase 
sequences  at  signals  1 and  2. 

i.  Express  the  inbound  and  outbound  bandwidths  as 
functions  of  the  offset  at  signal  2 using  the 
beginning  and  end  of  green  functions  obtained  in 
step  lb. 

ii.  Find  the  maximum  sum  of  the  inbound  and 
outbound  bandwidth  functions  obtained  in  step  2i. 

Step  3 . Select  the  phase  sequence  which  has  the  maximum  sum 
of  bands  as  determined  in  step  2 . 


CHAPTER  THREE 


EXPANDING  THE  ALGORITHM  TO  THE  GENERAL  ARTERIAL  HIGHWAY 


Optimal  phase  sequence  combinations  for  an  arterial 
highway  system  are  obtained  as  a natural  extension  of  the 
two-signal  solution  discussed  previously.  In  this  chapter 
an  algorithm  is  developed  to  obtained  the  maximum  sum  of 
bandwidths  by  successively  including  the  signals  of  the 
artery. 

For  an  arterial  highway  consisting  of  N signals, 

S, , S2, . . . , and  SN/  let  x;  denote  the  offset  at  signal  S;.  The 
discussion  in  chapter  2 above  shows  that  the  beginning  of 
the  green  at  each  signal  can  be  expressed  as  Xj  + c;  for 
some  constant  c*.  The  value  of  C;  depends  on  which  of  the 
four  phase  sequences  is  chosen  at  the  signal.  Likewise  the 
end  of  the  green  at  each  signal  can  be  expressed  as  X;  + k; 
for  some  constant  k;. 

Now,  let  TT(m,n)  represent  the  travel  time  from 
signal  m to  signal  n,  e;  = Xj  + kj  represent  the  end  of  the 
green  at  signal  S;,  and  b;  = x;  + c;  be  the  beginning  of 
the  green  at  signal  S;.  Then  E,  the  end  of  the  green  band 
on  the  artery  is  given  by 
E = min  {ei+TT(i,n)|  i = 1,  2,...  , N} 
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= min{xi+ci+TT(i,n)  | i = 1,  2,...,  N}  and  B,  the  beginning 
of  the  green  band  is  given  by 

B = max{bj+TT(i,n)  j i = 1,2,  . . . , N} 

= max  {Xi+kj+TT(i,n)  | i = 1,2,.  . .,  N} . 

These  values  are  used  to  see  that  the  bandwidth  on  the 
artery  is  expressible  as  a function  of  the  N variables, 
x,,  x2,  . . .,  xN.  Thus  it  has  been  demonstrated  that  the 
following  claim  is  true. 

CLAIM  1:  The  bandwidth  on  an  N signal  artery  is 

expressible  as 

R(x,,X2, . . . ,XN)  = min  {Xj+k;+TT  (i , n)  } 

- max  {Xi+Cj+TT  (i , n) 
where  x;  is  the  offset  at  signal  Si, 
k;  and  Cj  are  constants, 

TT ( i , N ) is  the  travel  time  from  signal  S-t  to 
signal  SN,  i = 1,2,3, ...,n. 

The  lemma  below  shows  that  a number  expressible  as 
a minimum  minus  a maximum  can  be  found  in  steps  by 
successively  adding  numbers  to  the  set  of  numbers  that  the 
minimum  and  maximum  are  taken  over.  Thus  Claim  1 above  and 
Lemma  1 below  will  show  that  the  bandwidth  on  an  artery  can 
be  found  by  first  finding  the  maximum  bandwidth  between  the 
first  and  last  signals,  then  successively  adding  constraints 
introduced  by  the  middle  signals  one  at  a time.  This  method 
of  computing  the  bandwidth  allows  the  computation  of  the 
maximum  sum  of  inbound  and  outbound  bandwidths  quite  simply 


37 


by  finding  the  maximum  value  of  a sum  of  two  piecewise 
linear  functions  as  illustrated  in  the  two-signal  case. 

LEMMA  1 • Let  A = { aj  , a2  , . . • , an } , B = ■{  bj , b2  / . . • , bn } , 

A'c  A,  and  B'c  B,  then 
min  A - max  B < min  A '-max  B'  . 
proof:  Since  min  A is  a minimum  over  more  numbers  than 
min  A'  it  follows  that  min  A < min  A' . Likewise 
max  B > max  B'  (since  the  maximum  is  taken 
over  more  numbers) . Thus  we  have  min  A < min  A' 
and  max  B > max  B'.  Multiplying  both  sides  of 
the  last  inequality  by  -1  gives  -max  B < -max 
B'.  Finally,  adding  the  two  inequalities  gives 
the  required  inequality, 

min  A - max  B < min  A'  - max  B'. 

COMPUTATION  OF  THE  GREEN  BAND  ON  AN  ARTERIAL  HIGHWAY 
Using  claim  1 and  lemma  1 above  the  bandwidth  on  an  artery 
can  be  calculated  by  computing  the  maximum  value  of  N-l 
piecewise  linear  functions  sequentially.  Denote 
min  { x,+k,+TT ( 1 , N)  ,xN+kN+TT(l,N)  } 

- max{x,+c1+TT(l,N)  , xN+cN+TT  ( 1 , N)  } 
by  R(x,,xn)  and  in  general  let  R(x,,x2,  . . . ,xi,xN)  denote 
min  {Xj+k^TTC^N)  , . . . , Xj+kj+TT ( i , N)  , xN+kN+TT ( 1 , N)  } 

- max{x1+c,+TT ( 1 , N)  , . . . , Xj+Cj+TT ( i , N)  , xN+cN+TT ( 1 , N)  } . 

Then  R(x,,xn)  is  the  maximum  bandwidth  between  signals  1 and 
N,  ignoring  the  remaining  signals.  R(x,,x2,  . . . ,Xj,xN)  is  the 
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bandwidth  on  the  artery  from  signals  1 to  i.  Lemma  1 shows 
that  R(x,  ,x2,  . • • (xN)  < R(x,,x2,  . • • /XN.2,xN} 

<.  R ( x, , x2 , . . . , xN_3 , xN ) 

< R(x, ,x2, . . . , x^jj , xN ) 

< ...  < ...<  R(x,,x2,xN)  < R(x,,xN). 

The  value  of  R(x,,xN)  can  be  computed  as  in  chapter  two  as  the 
maximum  of  a piecewise  linear  function  of  one  variable.  Let 
curmin,  = min{x,+k,+TT (1 , N)  , xN+kN+TT ( 1 , N)  } 
and  curmax,  = max{x,+c,+TT  ( 1 , N)  , xN+cN+TT  ( 1 , N)  } . 

Then  the  bandwidth  between  signals  1 and  N is  given  by 
R(x,  ,xN)  = curmin!  - curmax,.  Also, 

R(x,,x2,xn)=  min  {x,+k,+TT  ( 1 , N)  ,x2+k2+TT(2,N)  , xN+kN+TT  ( 1 , N)  } 

- maxlx^c^TTfl/N)  , x2+c2+TT  ( 1 , N)  } 

= min{x2+k2+TT ( 2, N)  , curmin,} 

- max  { x2+c2+TT  ( 2 , N ) , curmax, } . 

Similarly,  R(x,,x2, . . . ,x,,xN)  = min{xi+k,+TT(i,N)  , curmin,.,} 

- max { x,+c,+TT ( i , N ) , curmax,.,} 

where 

curmin,  = min{x,+k,+TT(l,N)  , ...  , x,+k,+TT  (i , N)  , xN+kN+TT(l,N)  } 
and 

curmax,  = max{x,+c,+TT ( 1 , N)  , . . . , x,+c,+TT(i,N)  , xN+cN+TT ( 1 , N)  } . 
After  N-l  iterations  of  this  process,  the  bandwidth  on  the 
artery  is  found  to  be 

R(x, , x2,  . . . , xN)  = min{xN.,+kN.,+TT (N-l , N)  , curminN.2} 

- max{xN.,+cN_,  +TT (N-l , N)  , curmaxN.2} 

= curminN_,  - curmaxN_,. 
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Thus  to  compute  the  bandwidth  on  the  artery,  first 
compute  curmin,  and  curmax, , then  compute  curminj  and 
curmaXj  for  i = 2 to  N-l  sequentially  as  maximum  values  of 
piecewise  linear  functions  as  in  chapter  2.  The  bandwidth 
is  then  curminN.,  - curmaxN.j . 

SUMMARY  OF  ALGORITHM  OPTIMIZING  AN  N SIGNAL  ARTERY 
Step  1.  Find  the  phase  sequence  combination (s)  and 
corresponding  offsets  which  yield  the  maximum  sum  of 
inbound  and  outbound  bands  between  signals  1 and  N. 
This  is  done  by  using  the  two  signal  arterial 
algorithm  presented  in  chapter  2. 

Do  the  following  for  each  phase  sequence  at  signal 
i,  i=  2 to  N-l. 

Step  2.  Fix  the  offset  and  phase  sequences  at  signals  1 and 
N.  Then  compute 

min{e,+TT(l,N)  , . . . , e;+TT(i,N)  , , eN} 

- max{b,+TT(l,N)  , . . . , bj+TT ( i , N)  , . . . , bN} 
for  each  of  the  four  phase  sequence  combinations. 
[Note  that  the  only  variable  in  the  above  expression  is  xi; 
the  offset  at  signal  SJ  . 

Step  3.  Select  the  maximum  of  the  four  values  computed  in 
step  2.  Fix  the  offset  , xf,  and  the  phase  sequence 
at  signal  i to  be  (one  of)  the  phase  sequences  which 
yields  the  maximum  sum  of  inbound  and  outbound 


bandwidths . 
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Example  Arterial  Highways 

To  see  how  the  algorithm  developed  in  this  chapter 
works,  two  example  arterial  highways  are  optimized  below. 
Appendix  B is  a listing  of  a computer  program  that  was 
written  and  used  as  an  aid  in  applying  the  algorithms  of 
this  research  to  the  examples. 

The  first  arterial  highway  to  be  optimized  consists  of 
three  signalized  intersections  along  an  east-west  corridor. 
The  distances  between  intersections  and  travel  times  are 
shown  in  table  3 . 1 and  the  green  requirement  per  movement  is 
shown  in  table  3.2. 


TABLE  3.1:  DISTANCES  BETWEEN  INTERSECTIONS  AND  TRAVEL  TIMES 

ON  EXAMPLE  1 


LINK 

LINK  LENGTH 

TRAVEL  TIME 

1 to  2 

1050  ft 

25  secs 

2 to  1 

1050  ft 

25  secs 

2 to  3 

975  ft 

23  sec 

3 to  2 

1000  ft 

2 0 sec 
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TABLE  3.2:  GREEN  TIME  REQUIREMENTS  BY  MOVEMENT  FOR  EXAMPLE 

1 


MOVEMENT 

SIGNAL  1 

SIGNAL  2 

SIGNAL  3 

Outbound  through 

30 

20 

35 

Outbound  left 

15 

25 

10 

Inbound  through 

25 

30 

25 

Inbound  left 

10 

15 

10 

Solution  to  Example  Arterial  Highway  1 

The  first  step  in  selecting  optimal  phase  sequence 
combinations  at  each  signal  of  the  arterial  highway  is  to 
determine  the  green  times  per  phase  for  each  of  the  four 
possible  left  turn  patterns.  These  assignments  are  shown 
below  in  table  3.3.  Using  the  green  time  per  phase,  the 
TABLE  3.3:  GREEN  TIME  ASSIGNMENTS  BY  MOVEMENT  FOR  EXAMPLE  1 


SIGNAL 

PHASE  SEQUENCE 

O 

I 

T 

L 

1 

Lead-LAG 

15 

10 

35 

0 

LAG-Lead 

15 

10 

35 

0 

LAG-LAG 

5 

0 

45 

0 

Lead-Lead 

5 

0 

45 

0 

2 

Lead-LAG 

25 

10 

15 

0 

LAG-Lead 

25 

10 

15 

0 

LAG -LAG 

10 

0 

25 

15 

Lead-Lead 

10 

0 

25 

15 

3 

Lead-LAG 

10 

10 

30 

0 

LAG-Lead 

10 

10 

30 

0 

LAG -LAG 

0 

0 

40 

10 

Lead-Lead 

0 

0 

40 

10 
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beginning  and  end  ofgreen  calculations  are  made  at  each 
intersection  for  each  phase  sequence. 

After  determining  the  green  time  per  phase,  the  next 
step  is  to  define  the  bandwidth  functions  and  finally  to 
calculate  the  maximum  sum  of  the  inbound  and  outbound  bands 
between  signals  1 and  3.  The  maximum  sum  for  each 
combination  of  phase  sequence  combinations  is  shown  in  table 
3.5  along  with  the  offsets  which  give  the  maximum  sum 


TABLE  3.4:  MAXIMUM  SUM  OF  BANDWIDTHS  BETWEEN  SIGNALS  1 AND 

3 FOR  EXAMPLE  1 


MAXIMUM  SUM 

OFFSET (S) 

PHASE  SEQUENCE 

OF  BANDS 

SIGNAL  2 

LAG-LAG/ LAG-LAG 

0 

0 

LAG-LAG/ Lead-Lead 

0 

0 

LAG-LAG/ Lead-LAG 

0 

0 

LAG-LAG/LAG-Lead 

73 

[48,55] 

Lead-Lead / LAG-LAG 

78 

[68,70] 

Lead-Lead/ Lead-Lead 

78 

[58,60] 

Lead-Lead/Lead-LAG 

80 

[60,65] 

Lead-Lead/LAG-Lead 

68 

[58,70] 

Lead-LAG  / LAG— LAG 

68 

[58,70] 

Lead-LAG/ Lead-Lead 

68 

[58,70] 

Lead— LAG  / Lead— LAG 

78 

[58,60] 

Lead— LAG / LAG— Lead 

58 

[48,70] 

LAG— Lead/ LAG— LAG 

80 

[58,60] 

LAG-Lead/Lead-Lead 

45 

[48,50] 

LAG-Lead/ Lead-LAG 

60 

[45,48] 

LAG-Lead/LAG-Lead 

60 

0 

Table  3 . 4 shows  that  two  phase  sequence  combinations 
yield  the  maximum  possible  sum  of  bands, lead-lead  /lead-lag, 
and  lag-lead/ lag-lag . Signal  2,  the  middle  signal,  is 
introduce  after  assigning  phase  sequence  lag-lead  and  offset 
0 at  signal  1 and  phase  sequence  lag-lag  and  offset  58  at 
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signal  3.  Table  3.5  shows  the  calculated  maximum  sum  of 
inbound  and  outbound  bands  for  each  of  the  four  phase 
sequence  combinations  at  signal  2.  Phase  sequence  lead-lag 
yields  a maximum  sum  of  bands  of  47.  Thus,  the  solution  to 
the  example  is  found  and  is  illustrated  in  table  3.6. 

Figure  3.1  shows  the  corresponding  time  space  diagram. 

TABLE  3.5:  MAXIMUM  SUM  OF  GREEN  BANDS  BY  PHASE  SEQUENCE  AT 

SIGNAL  2 OF  EXAMPLE  1 


PHASE  SEQUENCE 

OPTIMAL  OFFSET 

MAXIMUM  SUM  OF  BANDS 

Lag-Lag 

[78,93] 

25 

Lead-Lead 

[20,25] 

35 

Lag-Lead 

[40,53] 

47 

Lead-Lag 

[78,93] 

25 

TABLE  3.6:  OPTIMAL  SOLUTION  TO  EXAMPLE  1 


SIGNAL 

PHASE-  SEQUENCE 

OFFSET 

1 

Lag-Lead 

0 

2 

Lead-Lag 

40 

3 

Lag-Lag 

58 

INBOUND 

OUTBOUND 

BAND  12  SECS 
BAND  35  SECS 

Example  Arterial  Highway  #2 

The  next  arterial  example  is  what  might  be  considered  a 
worst  case  of  spacing  example.  This  arterial  highway  consists 
°f  five  signalized  intersections  along  an  east-west  corridor. 
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The  distances  between  intersections  and  travel  times  are 
shown  in  table  3.7  and  the  green  requirement  per  movement  is 
shown  in  table  3.8. 

TABLE  3.7:  DISTANCES  BETWEEN  INTERSECTIONS  AND  TRAVEL 

TRAVEL  TIMES  ON  EXAMPLE  2 


LINK 

LINK  LENGTH 

TRAVEL  TIME 

1 to  2 

1500  ft 

34 

secs 

2 to  1 

1500  ft 

34 

secs 

2 to  3 

660  ft 

15 

secs 

3 to  2 

660  ft 

15 

secs 

3 to  4 

660  ft 

15 

secs 

4 to  3 

660  ft 

15 

secs 

4 to  5 

1500  ft 

34 

secs 

5 to  4 

1500  ft 

34 

secs 

TABLE  3.8:  GREEN  TIME 

REQUIREMENTS 

2 

BY  MOVEMENT 

FOR  EXAMPLE 

MOVEMENT 

SIG 

1 SIG  2 

SIG 

3 

SIG 

4 SIG  5 

Outbound  through 

20 

20 

20 

20 

20 

Outbound  left 

10 

10 

10 

10 

10 

Inbound  through 

20 

20 

20 

20 

20 

Inbound  left 

10 

10 

10 

10 

10 

Solution  to  Example  Arterial  Highway  2 

First,  the  green  times  per  phase  for  each  of  the  four 
possible  left  turn  phase  sequences  are  computed.  These 
assignments  are  shown  below  in  table  3.9.  Using  the  green 
time  per  phase,  the  beginning  and  end  of  green  calculations 
are  made  at  each  intersection  for  each  phase  sequence. 
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TABLE  3.9: 

GREEN  TIME 

ASSIGNMENTS 
EXAMPLE  2 

BY  PHASE 

SEQUENCE 

FOR 

SIGNAL 

PHASE  SEQUENCE  0 

I 

T 

L 

1 

Lead-Lag 

10 

10 

20 

0 

Lag-Lead 

10 

10 

20 

0 

Lag-Lag 

0 

0 

30 

10 

Lead-Lead 

0 

0 

30 

10 

2 

Lead-Lag 

10 

10 

20 

0 

Lag-Lead 

10 

10 

20 

0 

Lag-Lag 

0 

0 

30 

10 

Lead-Lead 

0 

0 

30 

10 

3 

Lead-Lag 

10 

10 

20 

0 

Lag-Lead 

10 

10 

20 

0 

Lag-Lag 

0 

0 

30 

10 

Lead-Lead 

0 

0 

30 

10 

4 

Lead-Lag 

10 

10 

20 

0 

Lag-Lead 

10 

10 

20 

0 

Lag-Lag 

0 

0 

30 

10 

Lead-Lead 

0 

0 

30 

10 

5 

Lead-Lag 

10 

10 

20 

0 

Lag-Lead 

10 

10 

20 

0 

Lag-Lag 

0 

0 

30 

10 

Lead-Lead 

0 

0 

30 

10 

After  determining  the  green  time  per  phase,  the  next 
step  is  to  define  the  bandwidth  functions  and  finally  to 
calculate  the  maximum  sum  of  the  inbound  and  outbound  bands 
between  signals  1 and  3.  The  maximum  sum  for  each 
combination  of  phase  sequence  combinations  is  shown  in  table 
3.10  along  with  the  offsets  which  give  the  maximum  sum. 

Table  3.10  shows  that  phase  sequence  lead-lag/ lag-lead 
yields  the  maximum  possible  sum  of  bands  of  56.  Signal  2 is 
considered  next  after  assigning  phase  sequence  lead-lag  and 
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offset  0 at  signal  1 and  phase  sequence  lag-lead  and  offset 
28  at  signal  5. 

Table  3.11  shows  the  calculated  maximum  sum  of  inbound 
and  outbound  bands  for  each  of  the  four  phase  sequence 
combinations  at  signal  2.  Phase  sequence  lead-lead  yields  a 
maximum  sum  of  bands  of  56.  Next  signals  3 and  4 are 
introduced  into  the  optimization  process.  Eventually,  the 
solution  to  the  example  is  found  and  is  illustrated  in  table 
3.12. 

TABLE  3.10:  MAXIMUM  SUM  OF  BANDWIDTHS  BETWEEN  SIGNALS  1 AND 

5 FOR  EXAMPLE  2 


MAXIMUM  SUM 

OFFSET (S) 

Lag-Lag/ Lag-Lag 

44 

[22,38] 

Lag-Lag/ Lead-Lead 

44 

[12,20] 

Lag-Lag / Lead-Lag 

34 

[12,38] 

Lag-Lag/ Lag-Lead 

54 

[22,28] 

Lead-Lead/ Lag-Lag 

44 

[32,48] 

Lead-Lead/Lead-Lead 

44 

[22,38] 

Lead-Lead/ Lead-Lag 

34 

[22,48] 

Lead-Lead/ Lag-Lead 

54 

[32,38] 

Lead-Lag/ Lag-Lag 

54 

[32,38] 

Lead-Lag/ Lead-Lead 

54 

[32,38] 

Lead-Lag/ Lead-Lag 

44 

[22,38] 

Lead-Lag/ Lag-Lead 

56 

[28,32] 

Lag-Lead/ Lag-Lag 

34 

[22,48] 

Lag-Lead/ Lead-Lead 

30 

12 

Lag-Lead/ Lead-Lag 

30 

12 

Lag-Lead/ Lag-Lead 

44 

[22,38] 
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TABLE  3.11:  MAXIMUM  SUM  OF  GREEN  BANDS  BY  PHASE  SEQUENCE  AT 

SIGNAL  2 OF  EXAMPLE  2 


PHASE  SEQUENCE 

OPTIMAL  OFFSET 

MAXIMUM  SUM  OF  BANDS 

Lag-Lag 

34 

56 

Lead-Lead 

24 

56 

Lag-Lead 

[26,34] 

48 

Lead-Lag 

[24,32] 

48 

TABLE  3.12:  OPTIMAL  SOLUTION  TO  EXAMPLE  2 


SIGNAL 

PHASE-SEQUENCE 

OFFSET 

1 

Lead-Lag 

0 

2 

Lead-Lead 

40 

3 

Lag-Lead 

58 

4 

Lead-Lag 

4 

5 

Lag-Lead 

28 

INBOUND 

OUTBOUND 

BAND  8 secs 
BAND  30  secs 

48 


OUTBOUND  BANDWIDTH  - MINf  108, 103,  98}  - MAX)  58,  58,  63  } 
-98  - 63  - 35 


INBOUND  BANDW(TH  - MINff  145,143,115  ) - MAX)  103,  90, 100 } - 
-115-103-12 


FIGURE  3.1:  TIME  SPACE  DIAGRAM  OF  SOLUTION  TO  EXAMPLE  1 


CHAPTER  FOUR 


APPLICATIONS  TO  CLOSED  NETWORKS 

The  phase  sequence  optimization  problem  is  now  extended 
to  the  simple  closed  loop  consisting  of  four  arteries 
intersecting  at  four  signals. 


L 


s, 


Artery  1 


S2 


1 


FIGURE  4.1:  SIMPLE  CLOSED  LOOP 

Optimization  of  each  artery  of  the  loop  using  the 
procedures  of  chapter  2 above  yields  offsets  and  phase 
sequence  patterns  at  each  intersection  which  result  in  the 
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maximum  sum  of  bands  on  each  artery  of  the  loop.  The 
optimization  of  each  artery  yields  the  following  optimal 
assignments. 

Let  Xj  represent  the  offset  at  signal  Sif  EWPj  represent 
the  east-west  phase  pattern  at  signal  Su  and  NSP;  the  north 
south  phase  pattern  at  signal  Sj. 


Artery 

Phase  Pattern 

Offset 

Optimized 

Selected 

Selected 

Artery  1 

EWP(  and  EWP2 

x2 

Artery  2 

NSP2  and  NSP3 

x3 

Artery  3 

EWP3  and  EWP4 

x4 

Artery  4 

NSP,  and  NSP4 

X, 

The  values  selected  for  must  satisfy  the  following 
equations  and  inequalities. 

*2=  x,+  y2 , a2<  y2<  b2 

x3=  x2+  y3,  a3<  y3<  b3  Conditions  4.1 

x4=  x3+  y4,  a4<  y4<  b4 


xi=  x4+  y,  , a,<  y,<  b, 


[ Values  for  the  a/s  and  b/s  represent  the  endpoints  of  the 
interval  on  which  each  maximum  sum  is  attained] 

The  network  closure  constraint  requires  that  the 
offsets  are  selected  so  that  the  sum  of  the  offsets  around 
the  loop  should  be  an  integral  multiple  of  the  cycle  length. 
Thus  an  additional  requirement  on  the  x/s  is  that  Z Xj=mCL 

for  some  integer  m,  where  CL  represents  the  systemwide  cycle 
length. 
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LEMMA  2:  There  exists  x,  , i = 1,  2,  3,  4,  satisfying 

conditions  4.1  and  such  that  E x,=  mCL  if  and 

only 

if  Z as<  MCL  < Z b(. 

proof:  If  Z a,<  MCL  < Z b,  for  some  M,  then  consider 

case  1.  If  a2+  MCL  - Z a,  < b2 
Choose 

y2  = a2  + MCL  - Z a, 

y3  = a3 

y4  = a4 
y.=  a, 

then  aj<  yj<  bjfor  each  i, 


Xl  = [ 

2 MCL  + 

3a, 

+ a3  + 2a4]  / 4 

x2  = [ 

2MCL  - 

ai  - 

3 a3  2 a4  ] / 4 

x3  = [ 

2MCL  - 

a,+ 

a3-  2a4]  / 4 

x4  = [ 

2MCL  - 

a,  + 

a3+  2a4]  / 4 

then  Z x;  = 2*M*CL. 

Thus  if  m = 2M  then  Z x;=  mCL. 

case  2:  If  a2  + MCL  - Z a,>  b2,  then  choose 

y2=  b2and  let  k2=  a2+  MCL  - Z a;  -b2 
case  2A.  If  a3  + k2<  b3  then  choose 

y3=  a3+  k2 

= MCL  - a,  - a4  - b2 

y4=  a4 


Yi=  ai 
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and  x,  = 

: [2a,+  a4  - b2]/4 

x2  = 

[ 2a]  + a4=  3 b2 ] / 4 

x3  = 

■ [4MCL  - 2a,-  3a4-  b2]/4 

x4  = 

' [ 4 MCL  - 2a,  + a4  - b2]/4 

then  £ X;=  2MCL. 

case  2B . If  a3+  k2>  b3then  choose  y3=  b3  and 

let  k3=  b3  - a3-  k2  = b3  - MCL  + a,  + a4+  b2 
case  2Bi  . If  a4+  k3<  b4then  choose 


and  let 

y4=  a4+  k3 
y,=  a, 

x,  = [3MCL  - 4b2  - 3b3  - 2a4-  a,]/4 
x2=  [ 3MCL  - 3b3  - 2a4  - a,]/4 
X3=  [3MCL  + b3-  2a4  - a,]/4 
x4=  [MCL  + 5b3  + 6a4  + 3a,  + 4b2]/4 

then  £ x,=  2MCL 

case  2Bii.  If  a4  + k3  > b4then  choose 


and  let 

y4=  b4 

y,  = a,  + b4-b3-b2+  MCL  - a, -2a4 
x,  = [ 2MCL  - 3b2  - 2b3-b4]/4 
x2=  [ 2MCL  + b2-  2b3  - b4]/4 
x3=  [ 2MCL  + b2+  2b3  - b4]  /4 
x4=  [ 2MCL  +b2+  2b3  + 3b4]/4 

then  £ x,  = 2MCL. 
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Therefore  in  each  case  there  exist  X;  , i = 1, 2,3,4  which 
satisfy  the  conditions  required. 

Lemma  2 shows  that  a solution  which  is  optimal  on  all 
four  arteries  of  the  loop  exists  if  and  only  if  a multiple 
of  the  cycle  length  lies  between  the  sum  of  the  a/s  and  the 
sum  of  the  b/s.  The  following  example  illustrates  such  a 
case. 

EXAMPLE  3 A SIMPLE  CLOSED  LOOP 


1700 

ft 


1500  ft 


FIGURE  4.2:  ILLUSTRATION  OF  EXAMPLE  3 


TABLE  4.1:  GREEN  REQUIREMENT  BY  MOVEMENT  FOR  EXAMPLE  3 


SIGNAL 

NBL 

NBT 

SBL 

SBT 

EBL 

EBT 

WBL 

WBT 

1 

10 

15 

0 

20 

10 

20 

10 

20 

2 

0 

15 

0 

15 

10 

20 

10 

20 

3 

0 

15 

0 

15 

10 

20 

10 

20 

4 

10 

15 

0 

20 

10 

20 

10 

20 
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TABLE  4.2:  LINK  TRAVEL  TIMES  FOR  EXAMPLE  3 


LINK 

TT 

LINK 

TT 

1 TO  2 

15 

3 TO  4 

17 

2 TO  1 

15 

4 TO  3 

17 

2 TO  3 

34 

4 TO  1 

38 

3 TO  2 

34 

1 TO  4 

38 

The  algorithm  developed  in  chapter  two  yields  the 
maximum  sum  of  inbound  and  outbound  bands  on  the  four 
arteries  of  the  loop.  Table  4.3  below  shows  the  results  of 
applying  the  algorithm. 


TABLE  4.3:  MAXIMUM  SUM  OF  BANDWIDTHS  BY  ARTERY  OF  EXAMPLE  3 


MAXIMUM 

OPTIMAL  PHASE 

OFFSET 

ARTERY 

SUM  OF  BANDS 

PATTERN 

INTERVAL 

1 

40 

Lag-Lead/ Lag-Lag 

25  < y2  < 45 

2 

30 

Lag-Lead/ Lag-Lead 

y3  = 34 

3 

44 

Lag-Lead/ Lag-Lag 

27  < y4  < 43 

4 

24 

Lead-Lag/ Lag-Lead 

32  < y,  < 38 

Since  25  + 34  + 27  + 32  = 118  and  45  + 34  +43  + 38  = 

160,  it  follows  that  2*CL  = 120  lies  between  Z a;  and  Z b;. 
The  proof  of  claim  3 illustrates  how  the  offsets,  X;,  i = 
1,2,3,  4 can  be  chosen  to  satisfy  conditions  4.1.  Choose  y2  = 
25+120-118  = 27,  y3=  34,  y4  = 27,  and  y[  = 32.  Then  the  offsets 
X;,  at  each  signal  are  found  to  be  x,=  1,  x2=  28,  x3  = 2 , and  x4 
= 29  ( Z Xj  =60). 
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Thus  a solution  which  yields  the  maximum  sum  of  inbound 
and  outbound  bands  on  each  artery  of  the  network  is  found  and 
summarized  in  table  4.4  below. 

TABLE  4.4:  OPTIMAL  SOLUTION  TO  EXAMPLE  3 


SIGNAL 

EW  PHASING 

NS  PHASING 

OFFSET 

1 

Lag-Lead 

Lag-Lead 

1 

2 

Lag-Lag 

Lag-Lead 

28 

3 

Lag-Lead 

Lag-Lead 

2 

4 

Lag-Lag 

Lead-Lag 

29 

If  there  does  not  exist  M so  that  MCL  lies  between  2 aj 
and  2 b;then  at  least  one  artery  of  the  loop  will  attain  less 
than  its  absolute  maximal  sum  of  bands.  The  goal  in 
optimization  then  is  to  widen  the  intervals  in  a way  that  will 
include  a multiple  of  the  cycle  length  and  still  obtain  a sum 
of  bands  which  is  as  close  as  possible  to  the  maximum.  The 
appropriate  compromise  to  make  in  this  case  is  a matter  of 
judgement.  One  can  make  the  decision  based  on 

1.  Whether  one  direction  is  favored  over  others, 

2.  Whether  one  artery  is  less  important  than  others, 

3.  A desire  to  maintain  a minimum  sum  of  bands  throughout 

the  loop,  or 

4.  Some  other  engineering  decision. 

The  approach  taken  by  this  author  is  to  sacrifice  one 
artery  of  the  loop.  That  is,  if  the  engineer's  decision  is  to 
allow  one  artery  of  the  loop  to  attain  a suboptimal  sum  of 
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bands  in  order  to  satisfy  the  network  closure  constraint  (e.g. 
the  artery  with  the  lowest  maximum  sum  of  bands)  , then  a 
solution  yielding  the  maximum  sum  of  bands  on  the  remaining 
three  arteries  of  the  loop  is  possible.  Consider  the  set  of 
constraining  equations  4.1,  and  suppose  without  loss  of 
generality  that  artery  four  is  to  be  sacrificed.  Then  the  new 


set 

of 

constraints 

are 

x2  = 

Xi  + 

y2  / 

a2  < 

y2< 

b2 

*3  = 

x2  + 

Yj  / 

a3  < 

v| 

ro 

>i 

b3 

*4  = 

x3  + 

y4  / 

a4< 

y4< 

b4 

*1  = 

x4  + 

yi » 

v| 

o 

VI 

£ 

CL 

and 

s 

X:=  MCL 

for 

some  M. 

Constraints  4.2 


Claim  4:  There  exist  a multiple  of  CL  between 

a2  + a3  + a4  and  CL  + bj  + b3  + b4 
Proof:  Clearly 

a2  + a3  + a4  < CL  + b2  + b3  + b4  < 4 * CL . 

Case  1.  Now  suppose  4*CL  < a2+  a3+  a4  , then 
4 *CL  < a2+  a3+  a4<  CL  + b2+  b3  + b4<  4*  CL  . This 
implies  that  we  have  equality  throughout  the  inequality. 
That  is, 

a2+  a3  + a4=  CL  + b2+  b3  + b4=  4*  CL  . 

Thus  4*  CL  lies  in  the  interval. 

Case  2.  If  3*CL  < a2  + a3  + a4  then 
3*CL  < a,  + a,  + a4<  b2+  b,  + L,<  3*  CL  . 


•'4  — 


This  also  implies  equality  throughout  and  thus, 
a2  + a3  + a4  = CL  + b2  + b3  + b4  = 3 * CL  . 
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Thus  3*  CL  lies  in  the  interval. 

Case  3.  If  2*CL  < a2+  a3  + a4  then 

2 *CL  < a2  + a3+  a4<  b2+  b3  + b4<  3*  CL  . But  , if  3*CL 
is  between  a2+  a3  + a4and  b2+  b3+  b4then  case  2 above 
applies.  Meaning  that  3*CL  lies  in  the  interval. 
Otherwise  2*CL  lies  in  the  interval. 

Case  4.  If  CL  < a2  + a3  + a4,  then 
CL  < a2+  a3  + a4<  3*CL. 

In  this  case  either  2*CL  or  CL  lie  in  the  interval. 

Example  4.2  below  illustrates  how  the  proof  of  Claim  4.2 

can  be  used  to  "optimize"  a loop  where  the  absolute  maximum 

sum  of  bands  on  all  four  arteries  is  not  possible. 

EXAMPLE  4;  SIMPLE  CLOSED  LOOP  - SUBOPTIMAL  SOLUTION  ON  ONE 
ARTERY  REQUIRED 

TT (1,2)  = 17  SEC 


1700 

ft 


1500  ft 


750  ft 

TT (4,3)  = 21  SEC 


FIGURE  4.3:  ILLUSTRATION  OF  EXAMPLE  4 
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TABLE  4.5:  GREEN  REQUIREMENT  BY  MOVEMENT  FOR  EXAMPLE  4 


SIGNAL 

NBL 

NBT 

SBL 

SBT 

EBL 

EBT 

WBL 

WBT 

1 

10 

25 

25 

10 

10 

25 

25 

10 

2 

25 

10 

10 

25 

25 

10 

10 

25 

3 

10 

30 

10 

30 

10 

30 

10 

30 

4 

10 

30 

10 

30 

10 

30 

10 

30 

TABLE  4.6: 

LINK 

TRAVEL 

TIMES 

EXAMPLE  4 

LINK 

TT 

LINK 

TT 

1 TO 

2 

17 

3 TO 

4 

17 

2 TO 

1 

21 

4 TO 

3 

21 

2 TO 

3 

17 

4 TO 

1 

17 

3 TO 

2 

21 

1 TO 

4 

21 

The  algorithm  of  chapter  two  again  yields  the  optimum 
solution  on  each  artery  of  the  loop.  The  optimization  results 
are  shown  in  table  4.7  below. 


TABLE  4.7:  MAXIMUM  SUM  OF  BANDWIDTHS  BY  ARTERY  OF  EXAMPLE  4 


Maximum 

Optimal  Phase 

Sum  of  Bands 

Pattern 

offset  interval 

38 

lead- lag/ lead- lag 

27  < y,  < 29 

38 

lead-lag/ lead-lag 

27  < y2  < 29 

38 

lead-lag/ lead- lag 

27  < y3  < 29 

38 

lead-lag/ lead-lag 

27  < y4  < 29 
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Since  27  + 27  + 27  + 27  = 108  and  29  + 29  +29  + 29  = 
116,  it  is  clear  that  no  multiple  of  the  cycle  length  (60 
secs)  lies  between  108  and  116.  Thus  we  choose  one  artery  to 
be  "sacrificed."  Choosing  to  sacrifice  artery  4 the  offsets 
Xj,  at  each  signal  are  found  to  be  x,=  13,  x2=  42,  x3=  19,  and 
x4=  46  . (x,+  x2+  x3+  x4=  120). 

Thus  a solution  which  yields  the  maximum  sum  of  inbound 
and  outbound  bands  on  three  of  the  four  arteries  has  been 
found  and  is  summarized  in  the  table  4.8  below. 

TABLE  4.8:  OPTIMAL  SOLUTION  TO  EXAMPLE  4 


SIGNAL 

EW  PHASING 

NS  PHASING 

OFFSET 

SUM  OF 
BANDS 

1 

Lead-Lag 

Lead-Lag 

13 

38 

2 

Lead-Lag 

Lead-lag 

42 

23 

3 

Lead-Lag 

Lead-Lag 

19 

38 

4 

Lead-Lag 

Lead-Lag 

46 

23 

Example  4.3:  Daytona  Beach  Network 

To  show  how  the  procedure  described  above  can  be  used  in 
conjunction  with  existing  signal  timing  tools,  the  twelve 
signal  grid  network  used  in  the  TRANSYT-7F  user's  manual29  as 
an  example  will  now  be  optimized.  The  network  consists  of 
three  east  west  arterial  highways  and  four  north  south 
arterial  highways.  Figure  4.4  is  an  illustration  of  this 
network.  TRANSYT-7F  only  optimizes  cycle  length,  splits  and 
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offsets.  The  procedure  above  will  be  used  to  obtain  an  even 
better  solution  for  the  network  by  finding  phase  sequence 
patterns  at  each  intersection  which  using  the  given  splits, 
cycle  length  and  offsets  will  yield  the  greatest  sum  of 
inbound  and  outbound  bandwidths  on  the  arteries  of  the 
network.  Table  4.9  shows  each  optimal  value  as  determined  by 
TRANSYT-7F.  Using  these  values  the  procedure  above  yields 
optimal  phase  sequence  patterns  shown  in  table  4.10. 
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FIGURE  4.4:  ILLUSTRATION  OF  EXAMPLE  5 DAYTONA  BEACH  NETWORK 


61 


Table  4.9  TRANSYT-7F  OPTIMIZATION  RESULTS  ON  EXAMPLE  5 

DAYTONA  BEACH  NETWORK 


SIGNAL  SPLIT 

EW/NS 

OFFSET 

SIGNAL 

SPLIT 

EW/NS 

OFFSET 

1 

79/41 

0 

7 

34/86 

119 

2 

102/18 

107 

8 

68/52 

96 

3 

62/58 

42 

9 

66/54 

4 

4 

48/72 

41 

10 

80/40 

83 

5 

34/86 

115 

11 

78/42 

19 

6 

82/38 

98 

12 

82/38 

51 

TABLE  4.10:  OPTIMAL  PHASE  SEQUENCES  BY  SIGNAL  FOR  DAYTONA 

BEACH  NETWORK 


SIGNAL 

EAST-WEST 

SUBSEQUENCE 

NORTH-SOUTH 

SUBSEQUENCE 

1 

Lead-Lead 

Lag-Lead 

(ANY) 

2 

Lead-Lag 

Lag-Lead 

(ANY) 

3 

Lead-Lag 

Lag-Lead 

(ANY) 

4 

Lag-Lead 

Lag-Lag 

5 

Lag-Lead 

Lag-Lead 

(ANY) 

6 

Lag-Lead  (ANY) 

Lag-Lead 

(ANY) 

7 

Lag-Lead  (ANY) 

Lag-Lead 

(ANY) 

8 

Lag-Lead 

Lag-Lead 

(ANY) 

9 

Lag-Lead 

Lag-Lead 

(ANY) 

10 

Lag-Lead  (ANY) 

Lag-Lead 

(ANY) 

12 

Lag-Lead 

Lead-Lag 

ANY  - INDICATES  THAT  NO  PHASE  SEQUENCE  IS  BETTER  THAN  THE 
OTHERS 


Table  4.11  shows  a comparison  of  several  performance 
indexes  of  the  TRANSYT-7F  solution  and  the  solution  obtained 
using  the  procedure  of  this  research.  The  comparison  shows 
that  the  solution  obtained  using  the  new  procedure  improves 
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each  measure  of  performance  on  the  network  and  substantially 
improves  (10.62%  increase)  the  overall  sum  of  bandwidths. 


TABLE  4.11:  COMPARISON  OF  PERFORMANCE  MEASURES  ON  EXAMPLE  5 

DAYTONA  BEACH 


PERFORMANCE 

INDEX 

TRANSYT-7F  SOLUTION 
BANDS  NOT  OPTIMIZED 

SOLUTION  WITH 
BANDS  OPTIMIZED 

TRANSYT-7F  PI 

4113.21 

4100.22 

FUEL  CONSUMPTION 

9119.47 

9079.71 

SUM  OF  BANDWIDTHS  271.20 

300.00 

CHAPTER  FIVE 


CONCLUSIONS  AND  RECOMMENDATIONS 
Conclusions 

The  previous  chapters  have  shown  that  the  optimization 
of  the  signal  timing  components  in  networks  of  multiphase 
signals  can  be  accomplished  without  resulting  to  trial  and 
error  methods.  A closed  form  function  of  n variables  has 
been  defined  whose  dependent  variable  is  the  sum  of  inbound 
and  outbound  bandwidths  on  an  artery  and  whose  maximum  value 
is  the  maximum  sum  of  inbound  and  outbound  bandwidths  on  the 
arterial  highway.  The  functions  defined  can  be  used  to 
select  phase  sequence  combinations  which  yield  a maximal  sum 
of  inbound  and  outbound  bandwidths  on  the  arteries  of  a 
network  of  multiphase  signals. 

The  definition  of  the  bandwidth  function  requires  that 
the  cycle  length  and  green  allocation  at  each  intersection 
be  given.  The  values  are  then  used  to  define  the  bandwidth 
on  each  arterial  highway  as  a function  whose  independent 
variables  are  the  offsets  at  the  signals  of  the  network.  An 
analysis  of  the  function  showed  that  the  maximum  value  of 
the  function  can  be  found  by  evaluating  a few  simple 
algebraic  expressions.  The  maximum  sum  of  several  of  the 
bandwidth  functions  can  be  obtained  simply  and  quickly  (less 
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than  one  second  per  calculation)  on  a personal  computer. 
Using  these  maximum  values  a method  of  selecting  optimal 
phase  sequence  combinations  at  each  signal  of  a network  of 
multiphase  signals  has  been  demonstrated.  The  method 
selects  the  combination  of  phase  sequence  combinations  which 
give  the  maximum  sum  of  green  bands.  Thus,  the  procedure 
developed  is  a practical  tool  which  can  be  used  in 
conjunction  with  existing  programs  to  completely  optimize 
networks  of  multiphase  signals. 

The  algorithm  not  only  selects  optimal  phase  sequences 
for  arterial  highways  but  also  provides  a range  of  offsets 
associated  with  the  phase  sequences.  The  network  solution 
selects  optimal  phase  sequences  for  each  approach  of  the 
network . 

Drawing  on  the  fact  that  cycle  length,  splits,  and 
offset  optimization  algorithms  are  readily  available  to  the 
traffic  engineer  a final  step  in  the  network  optimization 
process  has  been  developed.  The  selection  procedure 
developed  has  the  advantage  of  requiring  very  little 
computer  resources  and  therefore  would  be  a cost  effective 
addition  to  the  network  optimization  process.  The  algorithm 
can  be  included  as  an  additional  step  in  network 
optimization  using  other  optimization  tools. 
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Recommendations 

Since  the  technique  developed  finds  all  phase  sequence 
combinations  which  yield  a maximum  sum  of  bands,  it  is 
suggested  that  further  research  be  done  on  ranking  the  four 
phase  sequence  combinations  and  on  selecting  one  sequence 
among  many  which  offer  the  same  or  nearly  the  same  sum  of 
bands . 

The  procedure  as  demonstrated  considers  all  four  of  the 
phase  sequence  as  viable  subsequences  at  each  signal  of  the 
network.  It  is  also  recommended  that  further  investigation 
on  the  possibility  of  excluding  one  or  more  of  the  phase 
sequences  at  some  signals  of  the  network  be  carried  out. 

The  fact  that  the  procedure  developed  can  be  used  to 
select  offsets  on  an  arterial  highway  suggests  that  further 
research  on  how  to  use  the  bandwidth  function  to  select 
offsets  in  a network  is  also  needed. 

A comparison  of  the  results  of  the  procedure  developed 
above  and  the  results  obtained  using  existing  optimization 
MAXBAND-8619  should  be  undertaken.  This  would  require 
significant  computer  resources  to  accommodate  MAXBAND-86  and 
that  the  procedure  be  more  efficiently  programmed  for 
computer  implementation. 

Finally,  it  is  recommended  that  the  algorithms 
developed  during  this  research  be  programmed  and  then  used 
in  conjunction  with  existing  programs,  such  as  TRANS YT, 
which  optimize  all  of  the  other  components  of  signal  timing. 
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If  this  algorithm  is  used  in  the  recommended  manner, 
networks  can  be  completely  optimized  using  less  computer 
time  and  thereby  spending  fewer  dollars. 


APPENDIX  A 


DEFINITIONS  OF  BANDWIDTH  FUNCTIONS 


For  each  of  the  sixteen  combinations  of  phase  sequences 
the  inbound  and  outbound  bandwidth  functions  are  of  the 
general  form 


Definitions  of  At,  Bj,  and  C;  for  each  of  the  sixteen  phase 
sequence  combinations  are  listed  below.  Definitions  are 
identified  using  the  phase  sequence  at  the  first  signal 
followed  by  the  phase  sequence  at  the  second  signal. 

The  following  symbols  are  used  in  defining  the  functions: 

Tj  = Number  of  seconds  of  green  allotted  to  the  through 
movement 

Lj  = Number  of  seconds  of  green  allotted  to  the  left  turn 
movement 

Oj  = Number  of  seconds  of  green  allotted  to  the  outbound 
through  and  outbound  left  turn  movement 

lj  = Number  of  seconds  of  green  allotted  to  the  inbound 
through  and  inbound  left  turn  movement  . 


1.  LAG -LAG /LAG -LAG  SEQUENCE  COMBINATION 

A,  = min  { TT  (1,2)  , TTtl^J+O^-Oj-Tj} 

B,  = max  { TT (1,2)  , TT(1, 2)  +0,+T1-02-T2} 
Ci  = T2+02-TT  ( 1 , 2 ) 


Ff(X) 


X + C; 

Aj  + Ci 

A:  + C:  + B:  - X 


X < Aj 
Aj  < X < Bj 
X > Bj. 
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1.  LAG -LAG /LAG -LAG  SEQUENCE  COMBINATION  CONTINUED 

A2  = min  { TT  (2,1)  , TT  (2,1)  +T1+I,-T2-I2} 

B2  = max  { TT  (2,1)  , TT  ( 2 , 1) +T,+I,-T2-I2} 

C2  = T2+I2-TT  (2,1) 

2.  LAG-LAG/LEAD-LEAD  SEQUENCE  COMBINATION 

A,  = min  (TT(1,2) -I2-L2,  TT ( 1 , 2 ) -I2-L2+0,+T1-02-T2} 

B,  = max  {TT  (1,2)  -I2-L2 , TT  (1,2)  -I2-L2+0,+T,-02-T2} 

C,  = Ij+I^+O^^-TT  ( 1 , 2 ) 

A2  = min  {TT(2,l)-L2-02/  TT  (2,1)  -L2-02+T,+I1-T2-I2} 
B2  = max  {TT(2, 1) -Lj-Oj,  TT  (2 , 1)  -L2-02+T,+I,-T2-I2} 
C2  = I2+L2+02+T2-TT  ( 2 , 1 ) 

3.  LAG-LAG/ LEAD-LAG  SEQUENCE  COMBINATION 

a . Outbound 

A,  = min  {TT(1,2),  TT (1,2)  +0,+T,-02-T2} 

B,  = max  { TT  (1,2)  , TT  ( 1 , 2 ) +01+T1-02-T2} 

C,  = 02+T2-TT  ( 1 , 2 ) 

b.  Inbound 

A2  = min  {TT  (2,1)  -02 , TT  (2,1)  -02+T,+I1-T2-I2 
B2  = max  { TT (2,1)  -02 , TT(2 , 1) -02+T1+I,-T2-I2 
C2  = 02+I2+T2-TT  (2,1) 
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4.  LAG -LAG /LAG -LEAD  SEQUENCE  COMBINATION 

a.  Outbound 

A,  = min  {TT(1,2)-I2,  TT  ( 1 , 2 ) -I2+T,+0,-T2-02} 

B!  = max  {TT  ( 1 , 2 ) -I2 , TT  (1,2)  -I2+T1+0,-T2-02} 

C,  = I2+02+T  (2,4)  -TT  (1,2) 

b.  Inbound 

A2  = min  {TT (2,1)  , TT (2,1)  +T,+I,-T2-I2} 

B2  = max  {TT (2,1)  , TT(2 , 1)  +T,+I1-T2-I2} 

C2  = I2+T2-TT  (2,1) 

5.  LEAD-LEAD/ LAG-LAG  SEQUENCE  COMBINATION 

a.  Outbound  bandwidth 

A,  = min  { TT (1,2)  +L,+I, , TT (1,2)  +L1+I,+01+T1-02-T2} 

B,  = max  { TT(1,2)+L,+I,/  TT  (1,2)  +L,+I,+01+T1-02-T2} 

C,  = T2+02-L,-l!-TT  (1,2) 

b.  Inbound  bandwidth 

A2  = min  {TT  (2,1)  +L1+0I , TT(2 , 1) +L,+01+I,+T,-I2-T2  } 

B2  = max  {TT(2,  l)+L,+0!,  TT(2 , 1) +L1+01+I,+T,-I2-T2  } 

C,  = T(2,l) +1(2,1) -L (1,2) -O (1,2) -TT (2,1) 

6.  LEAD -LEAD /LEAD -LEAD  SEQUENCE  COMBINATION 

a . Outbound 

A,  = min  {TT  (1,2)  +L,+I1-L2-I2 , TT  (1,2)  +L1+I1-L2-I2+0,+T1-02-T2 

B,  = max  {TT  (1,2)  +L,+I1-L2-I2 , TT  (1,2)  +L,+I1-L2-I2+0,+T1-02-T2 

C,  L2+02+I2+T2-L,-I,-TT  (2,1) 
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b.  Inbound 

A2  = min{TT(2, 1)  +L,+0,-L2-02,  TT  (2 , 1) +L,+01-L2-02+I1+T,-I2-T2} 
B2  = max{TT (2,1)  +L,+0,-L2-02,  TT (2,1)  +L1+0,-L2-02+I,+T1-I2-T2} 
Cj  = L2+02+I2+T2-L,-01-TT  (2,1) 

7.  LEAD-LEAD /LEAD-LAG  SEQUENCE  COMBINATION 

a . Outbound 

A,  = min  {TT (1,2)  +L,+I1  , TT (1,2)  +L1+I1+0,+T,-02-T2} 

Bt  = max  {TT  (1,2)  +L,+Ij  , TT  ( 1 , 2 ) +L,+I1+01+T,-02-T2} 

Cj  = 02+T2-L,-I,-TT(l,2) 
b.  Inbound 

A2  = min  {TT(2,l)+L1+0,-02,  TT  (2 , 1) +L1+01-02+I1+T,-I2-T2} 

B2  = max  {TT  (2,1)  +L,+0,-02;  TT  (2 , 1) +L,+0,-02+I1+T1-I2-T2} 

Cj  = 02+I2+T2-L,-0,-TT  (2,1) 

8.  LEAD-LEAD/LAG-LEAD  SEQUENCE  COMBINATION 

a.  Outbound 

A,  = min  {TT(1,2)+L1+I1-I2,TT(1,2)+L1+Ii-I2+0,+T1-02-T2} 

B,  = max  {TT(1,2)+L1+I1-I2,TT(1,2)+L1+Ii-I2+0,+T,-02-T2} 

C,  = I2+02+T2-L,-I,-TT  (1,2) 
b.  Inbound 

A2  = min  {TT  (2,1)  +L,+0, : , TT  (2,1)  +L1+0,+I,+T,-I2-T2} 

B2  = max  {TT(2f  l)+L,+0,:  ,TT(2 , 1) +L,+01+I,+T1-I2-T2} 

C2  = I2+T2-L,-01-TT(2,1) 
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9.  LEAD -LAG /LAG -LAG  SEQUENCE  COMBINATION 

a.  Outbound 

A,  = min  { TT  (1,2)  , TT(1,2)+0,+Ti-02-T2} 

B,  = max  {TT  (1,2)  , TT  ( 1 , 2 ) +0,+T,-02-T2} 

Cj  = 02+T2-TT  ( 1 , 2 ) 

b.  Inbound 

A2  = min  {TT (2,1)  +0, , TT (2,1)  +01+T,+I1-T2-I2} 

B2  = max  {TT  (2,1)  +0, , TT  ( 2 , 1)  +0,+T1+I1-T2-I2} 

C2  = T2+I2-0,-TT  (2,1) 

10.  LEAD-LAG/LEAD-LEAD  SEQUENCE  COMBINATION 

a.  Outbound 

Ai  = min  {TT  ( 1 , 2 ) -I2,  TT  (1,2)  -I2+0,+T,-02-T2} 

B,  = max  {TT  ( 1 , 2 ) -I2 , TT  (1,2)  -Ij+Oj+Tj-Oj-Tj} 

C,  = I2+02+T2-TT  ( 1 , 2 ) 

b.  Inbound 

A2  = min  { TT  (2,1)  +0, -02 , TT  (2,1)  +01-02+T1+I,-T2-I2} 
B2  = max  {TT  (2,1)  +0,-02,  TT  (2,1)  +0,-02+T,+I1-T2-I2} 
C2  = 02+T2+I2-0,-TT  (2,1) 

11.  LEAD-LAG/ LEAD-LAG  SEQUENCE  COMBINATION 
a.  Outbound 

Aj  = min  { TT  (1,2)  , TT  (1,2)  +0,+T,-02-T2} 

B,  = max  {TT  (1,2)  , TT  (1,2)  +01+T,-02-T2} 

C,  = 02+T2-TT  ( 1 , 2 ) 
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b.  Inbound 

A2  = min  {0,-02+TT (2,1)  , T1+I,-I2-T2+01-02+TT(2,l)  } 
B2  = max  {0,-02+TT  (2,1)  , T^^-Ij-T^O^+TT  (2 , 1)  } 
C2  = 02+I2+T2-0,-TT  (2,1) 

12.  LEAD-LAG /LAG-LEAD  SEQUENCE  COMBINATION 

a.  Outbound 

A,  = min  {TT  (1,2)  -I2 , TT  (1,2)  -I2+0,+T1-T2-02} 

B,  = max  {TT  ( 1 , 2 ) -I2 , TT  (1,2)  -I2+0,+T,-T2-02} 

Cj  = 02+T2+I2-TT  ( 1 , 2 ) 

b.  Inbound 

A2  = min  {TT(2,l)+0,,  TT (1,2)  +0,+T,+I1-T2-I2} 

B2  = max  {TT  (2,1)  +0, , TT  ( 1 , 2 ) +0,+T,+I1-T2-I2} 

C2  = T2+I2-0,-TT  (2,1) 

13.  LAG-LEAD /LAG-LAG  SEQUENCE  COMBINATION 

a.  Outbound 

A,  = min  {TT  (1 , 2 ) +1, , TT  ( 1 , 2 ) +I1+0i+T1-02-T2} 

B, =  max  {TT  (1,2)  +1, , TT  (1,2)  +I1+01+T1-02-T2} 

C,  = T2+02-I,-TT  (1,2) 

b.  Inbound 

A2  = min  { TT  (2,1)  , TT(2 , 1)  +T,-I2-T2> 

B2  = max  { TT  (2,1)  , TT  (2 , 1) +T,-I2-T2} 

C2  = T2+I2-TT  (2,1) 
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14.  LAG- 

-LEAD /LEAD-LEAD  SEQUENCE  COMBINATION 

a. 

Outbound 

A,  = 

min  {TT  (1,2)  +I1-L2-I2,  TT  ( 1 , 2 ) +I1-L2-I2+01+T,-02-T2} 

B,  = 

max  {TT(1,2)+I,-L2-I2/  TT (1,2)  +I1-L2-I2+01+T,-02-T2} 

c,  = 

L2+I2+02+T2-I1-T,-TT  (1,2) 

b. 

Inbound 

ii 

<T 

min  {TT(2,l)-L2-02,  TT(2,l)-L2-02+I,+T,-I2-T2} 

b2  = 

max  {TT(2,l)-L2-02,  TT  (2 , 1) -L2-02+I,+T1-I2-T2} 

c2  = 

L2+02+I2+T2-TT  (2,1) 

15.  LAG-LEAD /LEAD-LAG  SEQUENCE  COMBINATION 


a. 

Outbound 

A,  = 

min  {TT  (1,2)  -I, , TT  (1 , 2)  -Ij+Oj+T^-TJ 

B,  = 

max  {TT  (1,2)  -Ij , TT  ( 1 , 2 ) -I^O^T,^-^} 

c,  = 

= 02+T2-I,-TT  (1,2) 

b. 

Inbound 

ii 

< 

= min  { TT  (2,1)  -02 , TT  (2 , 1) -02+I,+T1-I2-T2} 

b2  = 

= max  { TT  (2,1)  -02 , TT  (2 , 1) -02+I,+T,-I2-T2> 

c2  = 

= 02+I2+T2-TT  (2,1) 

16.  LAG-LEAD /LAG-LEAD  SEQUENCE  COMBINATION 


a. 

Outbound 

A,  = 

min  {TT  (1,2)  +I,-I2 , TT  (1,2)  +I1-I2+01+T,-02-T2} 

B,  = 

max  {TT(1,2)+I1-I2,  TT (1,2)  +I,-I2+01+T1-02-T2} 

c,  = 

I2+02+T2-I,-TT  (1,2) 

b.  Inbound 
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A2  = min  { TT  (2,1)  , TT  ( 2 , 1) +I,+T,-I2-T2} 
B2  = max  { TT (2,1)  , TT(2 , 1) +I,+T,-I2-T2} 
C2  = I2+T2-TT  (2,1) 


APPENDIX  B 
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315 
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PROGRAM  LISTING 

INPUT  "NUMBER  OF  SIGNALS  ";N 
INPUT  "CYCLE  LENGTH  ";C 
REM 

*************************************************** 
REM  * 1 = LAG-LAG,  2 = LEAD-LEAD,  3 = LEAD-LAG,  AND 
4 = LAG-LEAD* 

REM 

************************************************** 
DIM  OBL (N) : DIM  TT (N, N) : DIM  IBL(N):DIM  OBT(N):DIM 
IBT(N) : DIM  CL (N) 

DIM  A(2 , 16) :DIM  B(2,16):DIM  C(2,16):DIM 
X (2 , N, 4 ) : DIM  Y(2,N,4):DIM  Z(2,N,4) 

DIM  OFFSET (N) : DIM  OPTPAT (N) : DIM  FL(N,4):DIM 
PS ( 16 ) : DIM  PZ (N , 4 ) 

DIM  V(2 , 16) : DIM  U(2,16):DIM  W(2,N,4):DIM 
Q ( 2 , N , 4 ) : DIM  FA(16) 

DIM  MAXB(16) ,OFFB(16) ,OFFE(16) :DIM  VARB (N , 4 ) : DIM 
VARE(N,4) :DIM  GRBD (N, 4 ) 

DIM  BOB (N) : DIM  EOB(N):DIM  BIB(N):DIM  EIB(N) 

DIM 

0 (N , 4 ) ,I(N,4) , T (N , 4 ) ,L(N,4) , OB (N , 4 ) , OE (N,4) ,IB(N,4) 
, IE (N , 4 ) 

ICURMIN=10000 : OCURMIN=10000 : ICURMAX=0 : OCURMAX=0 

FOR  Q = 1 TO  N 

FOR  J=1  TO  N 

TT ( J, Q)  = 0 

NEXT  J 

NEXT  Q 

FOR  1=1  TO  N-l 

PRINT  "TRAVEL  TIME  SIGNAL" ; I ; "TO  SIGNAL" ;I+1; 

INPUT  TT( 1,1+1) 

PRINT  "TRAVEL  TIME  SIGNAL" ; 1+1 ; " TO  SIGNAL" ;I; 

INPUT  TT( 1+1,1) 

NEXT  I 

FOR  Z=1  TO  N-2 

FOR  A=  Z+2  TO  N 

FOR  B = Z TO  A-l 

TT ( Z , A) =TT ( Z , A) +TT (B , B+l ) 

TT (A, Z) =TT (A, Z) +TT (B+l , B) 

NEXT  B 
NEXT  A 
NEXT  Z 
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FOR  J=1  TO  N 
FOR  Z = 1 TO  N 

IF  TT ( Z , J) >C  THEN  TT ( Z , J) =TT ( Z , J) -C : GOTO  390 

IF  TT ( J , Z) >C  THEN  TT ( J, Z) =TT ( J, Z) -C: GOTO  400 

NEXT  Z 

NEXT  J 

FOR  1=1  TO  N 

PRINT  "OUTBOUND  LEFT, INBOUND  LEFT , OUTBOUND 
THRU, INBOUND  THRU  FOR  SIGNAL"; I 
INPUT  OBL(I),  IBL(I) , OBT (I) , IBT (I) 

INPUT  "NUMBER  SECONDS  GREEN  ALLOCATED  TO 
ARTERIAL  AT  THIS  SIGNAL  ";CL(I) 

NEXT  I 

FOR  1=1  TO  N 

0(1,3) =OBL ( I ) 

1(1,3) =IBL (I) 

T ( 1 , 3 ) =CL ( I ) -O (1,3) -I  (1,3) 

IF  T (I , 3 ) <0  THEN  T (I , 3 ) =CL (I) :0 (I , 3 ) =0 : I (I , 3 ) =0 
0(1,4) =0 (1,3) 

1(1,4) =1 (1,3) 

T ( I , 4 ) =T (1,3) 

IF  OBL(I) <IBL(I)  THEN  L ( I , 2 ) =OBL ( I ) ELSE 
L ( I , 2 ) =IBL ( I ) 

IF  IBL ( I ) -L ( 1 , 2 ) >0  THEN  I ( 1 , 2 ) = ( IBL ( I ) -L ( 1 , 2 ) ) 
ELSE  1(1,2) =0 

IF  OBL(I) -L(I,2) >0  THEN  0 ( I , 2 ) = (OBL ( I) -L ( I , 2 ) ) 
ELSE  0 ( I , 2 ) =0 

T ( I , 2 ) =CL (I) -L (I,2)-I(I,2) -O (1,2) 

IF  T ( I , 2 ) <0  THEN 

T ( 1 , 2 ) =CL ( I ) : 0 (1 , 2 ) =0 : 1 (1 , 2 ) =0 : L (1 , 2 ) =0 
L(I,1)=L(I,2) : I (1,1) =1(1,2) : 0(1,1) =0(1,2 ) :T (1,1 
) =T (1,2) 

NEXT  I 

REM  LELE/LALA  ROUTINE 

IF  0(1,2) +T ( 1 , 2 ) >0 ( N , 1 ) +T ( N , 1 ) THEN 

A ( 1 , 5 ) =TT ( 1 , N) +L ( 1 , 2 ) +1 ( 1 , 2 ) :B(1,5) =TT (1,N)+L(1,2)+ 
1(1,2) +0 (1,2) +T (1,2) -O ( N , 1 ) -T ( N , 1 ) ELSE 
A ( 1 , 5 ) =TT ( 1 , N ) +L ( 1 , 2 ) +1 ( 1 , 2 ) +0 (1,2) +T (1,2) -O ( N , 1 ) -T 
(N,  1)  : B ( 1, 5)  =TT  (1 , N)  +L(  1 ,2) +1(1, 2) 

C ( 1 , 5 ) =T (N, 1) +0 (N, 1) -L (1,2)-I(1,2) -TT ( 1 , N) 

IF  1(1,2) +T (1,2)>I(N,1) +T ( N , 1 ) THEN 

A ( 2 , 5 ) =C-TT ( N , 1 ) +L (1,2) +0 (1,2) : B (2 , 5) =C-TT (N, 1) +L(1 
, 2 ) +0 (1,2) +1 (1,2) +T (1,2)-I(N,1) -T (N , 1)  ELSE 
A ( 2 , 5 ) =C-TT (N , 1 ) +L (1,2) +0 (1,2) +1 (1,2) +T (1,2)-I(N,1) 
-T  (N,  1)  : B (2 , 5)  =C-TT  (N,  1)  +L (1,2)  +0  (1,2) 

C ( 2 , 5 ) =T (N,1)+I(N,1) -L (1,2) -O (1,2) +TT (N , 1) -C 

REM  LELE/LELE  ROUTINE 

IF  0(1,2) +T (1,2) >0 ( N , 2 ) +T ( N , 2 ) THEN 

A ( 1 , 6 ) =TT ( 1 , N) +L ( 1 , 2 ) +1 ( 1 , 2 ) -L (N , 2 ) -I (N , 2 ) : 

B ( 1 , 6 ) =A (1,6) +0 (1,2) +T (1, 2) -O ( N , 2 ) -T ( N , 2 ) ELSE 
B ( 1 , 6 ) =TT ( 1 , N) +L (1,2)+I(1,2) -L (N , 2 ) -I (N, 2 ) :A(1,6)=B 
(1,6)  +0  (1,2)  +T  (1,2)  -O  (N,  2 ) -T  (N,  2 ) 
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760 

770 


780 

790 

800 


810 

820 


830 

840 

850 


860 

870 


880 

890 

900 


910 

920 


930 

940 

950 


C ( 1 , 6 ) =L  (N , 2 ) +0 (N,2)+I(N,2) +T (N , 2 ) -L (1,2)-I(1,2) -TT 
(1,N) 

IF  1(1,2) +T (1,2)>I(N,2) +T ( N , 2 ) THEN 
A ( 2 , 6 ) =C-TT (N , 1 ) +L (1,2) +0 (1,2) -L (N, 2 ) -O  (N , 2 ) : 

B ( 2 , 6 ) =A (2,6) +1 (1,2) +T (1,2)-I(N,2) -T (N , 2 ) ELSE 
B ( 2 , 6 ) =C-TT  ( N , 1)  +L  (1,2)  +0  (1,2)  -L (N , 2 ) -O  (N,2)  : A ( 2 , 6 ) 
=B (2,6)+I(l,2) +T (1,2)-I(N,2) -T (N , 2 ) 

C ( 2 , 6 ) =L  (N , 2 ) +0  (N,  2 ) +1  (N , 2 ) +T  (N , 2 ) -L  (1,2)  -O  (1,2)  +TT 
(N,l)-C 

REM  LELE/LELA  ROUTINE 

IF  0(1,2) +T (1,2) >0 ( N , 3 ) +T ( N , 3 ) THEN 

A ( 1 , 7 ) =TT ( 1 , N) +L (1,2) +1(1, 2) : B ( 1 , 7 ) =A (1,7) +0(1,2)+ 

T ( 1 , 2 ) -O (N , 3 ) -T (N, 3 ) ELSE 

B ( 1 , 7 ) =TT ( 1 , N ) +L ( 1 , 2 ) +1 ( 1 , 2 ) : A ( 1 , 7 ) =B (1,7) +0 (1,2) +T 

(1.2)  -O (N , 3 ) -T (N , 3 ) 

C ( 1 , 7 ) =0 (N , 3 ) +T (N , 3 ) -L (1,2)-I(1,2) -TT ( 1 , N) 

IF  I ( 1 , 2 ) +T ( 1 , 2 ) >1 (N , 3 ) +T (N, 3 ) THEN 

A ( 2 , 7 ) =C-TT ( N , 1 ) +L (1,2) +0 (1,2) -O (N,3) :B(2,7) =A (2,7) 
+1 ( 1 , 2 ) +T ( 1 , 2 ) -I (N , 3 ) -T (N, 3 ) ELSE 

B ( 2 , 7 ) =C-TT (N , 1) +L (1,2) +0 (1,2) -O (N , 3 ) :A(2,7)=B(2,7) 
+1 (1,2) +T ( 1 , 2 ) -I (N, 3 ) -T (N , 3 ) 

C ( 2 , 7 ) =0  (N , 3 ) +1  (N , 3 ) +T  (N , 3 ) -L  (1,2)  -O  (1,2)  +TT  (N,  1)  -C 

REM  LELE/LALE  ROUTINE 

IF  0(1,2) +T (1,2) >0 ( N , 4 ) +T ( N , 4 ) THEN 

A ( 1 , 8 ) =TT ( 1 , N) +L (1,2) +1 (1,2) -I (N, 4 ) :B(1,8)=A(1,8)+ 

0(1,2) +T (1,2) -O (N , 4 ) -T (N , 4 ) ELSE 

B ( 1 , 8 ) =TT ( 1 , N) +L (1,2) +1 (1,2) -I  (N , 4 ) : A ( 1 , 8 ) =B ( 1 , 8 ) +0 

(1.2)  +T (1,2) -O ( N , 4 ) -T ( N , 4 ) 

C ( 1 , 8 ) =1  (N , 4 ) +0  (N , 4 ) +T  (N , 4 ) -L  (1,2)  -I  (1,2)  -TT  ( 1 , N) 

IF  1(1,2) +T (1,2)>I(N,4) +T ( N , 4 ) THEN 

A ( 2 , 8 ) =C-TT (N, 1) +L (1,2) +0  (1,2)  : B ( 2 , 8 ) =A (2, 8) +1(1, 2) 
+T(1,2)-I(N,4)-T(N,4)  ELSE 

B ( 2 , 8 ) =C-TT ( N , 1 ) +L (1,2) +0 (1,2) : A ( 2 , 8 ) =B ( 2 , 8 ) +1 ( 1 , 2 ) 
+T(1,2)-I(N,4)-T(N,4) 

C ( 2 , 8 ) =1  (N , 4 ) +T  (N , 4 ) -L  (1,2)  -O  (1,2)  +TT  (N , 1)  -C 

REM  LALA/LALA  ROUTINE 

IF  T(l, 1) +0(1, 1) >T(N, 1) +0(N, 1) 

A ( 1 , 1 ) =TT ( 1 , N ) : B (1,1) =TT(1,N) +0(1,1) +T( 1,1) -0(N,1)- 
T (N, 1)  ELSE 

A ( 1 , 1 ) =TT ( 1 , N ) +0 (1,1) +T(1, 1) -0(N, 1) -T(N, 1) :B(1, 1)=T 
T ( 1 , N) 

C(l, 1)=T(N, 1) +0(N, 1) -TT ( 1 , N) 

IF  T(l,l)+I(l,l) >T (N, 1) +1 (N, 1)  THEN 

A ( 2 , 1) =C-TT (N , 1) : B (2 , 1) =C-TT (N, 1) +T ( 1 , 1) +1 ( 1 , 1) -T (N 
, 1) -I (N, 1)  ELSE 

A ( 2 , 1 ) =C-TT (N , 1 ) +T (1,1)+I(1,1) -T(N,1)-I(N,1) :B(2,1) 
=C-TT(N, 1) 

C (2 , 1) =T (N, 1) +1 (N, 1) +TT (N, 1) -C 

REM  LALA/LELE  ROUTINE 

IF  T ( 1 , 1 ) +0 (1,1) >T ( N , 2 ) +0 ( N , 2 ) THEN 

A ( 1 , 2 ) =TT (1,N)-I(N,2) -L (N , 2 ) : B ( 1 , 2) =TT (1,N) -I (N, 2) - 

L ( N , 2 ) +0 (1,1) +T (1,1) -0 ( N , 2 ) -T ( N , 2 ) ELSE 
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A ( 1 , 2 ) =TT (1,N)-I(N,2) -L (N , 2 ) +0 (1,1) +T (1,1) ~0 (N, 2 ) -T 
(Nf2) : B ( 1 , 2 ) =TT (1,N)-I(N,2) -L (N , 2 ) 

960  C ( 1 , 2 ) =1 ( N , 2 ) +L ( N , 2 ) +0 ( N , 2 ) +T ( N , 2 ) -TT ( 1 , N ) 

970  IF  T(1,1)+I(1,1)>T(N,2)+I(N,2)  THEN 

A ( 2 , 2 ) =C-TT (N , 1 ) -L (N , 2 ) -0 (N,  2 ) : B (2 , 2) =C-TT (N, 1) -L(N 
,2)-0(N,2)+T(l,l)+I(l,l)-T(N,2)-I(N,2)  ELSE 
A(2 , 2) =C-TT (N,  1) -L(N,2) -0(N,2)+T(1, 1) +1(1,1) -T(N,2) 
-I (N, 2 ) : B ( 2 , 2 ) =C-TT  (N , 1)  -L (N, 2 ) -0 (N, 2 ) 

980  C ( 2 , 2 ) =1 (N , 2 ) +L (N , 2 ) +0 (N , 2 ) +T (N, 2 ) +TT (N , 1)-C 

990  REM  LALA/LELA  ROUTINE 

1000  IF  T ( 1 , 1 ) +0 (1,1) >T ( N , 3 ) +0 ( N , 3 ) THEN 

A ( 1 , 3 ) =TT (1,N):B(1,3) =TT (1,N)+0(1,1)+T(1,1)-0(N,3)- 
T (N , 3 ) ELSE 

A ( 1 , 3 ) =TT ( 1 , N ) +T (1,1) +0 (1,1) -T ( N , 3 ) -O ( N , 3 ) : B ( 1 , 3 ) =T 
T(1,N) 

1010  C ( 1 , 3 ) =0 ( N , 3 ) +T ( N , 3 ) -TT ( 1 , N ) 

1020  IF  T ( 1 , 1 ) +1 (1,1) >T (N , 3 ) +1 (N, 3 ) THEN 

A ( 2 , 3 ) =C-TT ( N , 1 ) -O ( N , 3 ) : B (2 , 3 ) =C-TT  (N,  1) -O  (N,  3 ) +T  (1 
,1)+I(1,1)-T(N,3)-I(N,3)  ELSE 

A(2,3)=C-TT(N,1)-0(N,3)+T(1,1)+I(1,1)-T(N,3)-I(N,3) 

: B (2 , 3 ) =C-TT (N, 1) -0(N,3) 

1030  C ( 2 , 3 ) =0 (N , 3 ) +1 (N , 3 ) +T (N , 3 ) +TT (N , 1) -C 

1040  REM  LALA/LALE  ROUTINE 

1050  IF  T ( 1 , 1 ) +0 (1,1) >T ( N , 4 ) +0 ( N , 4 ) THEN 

A ( 1 , 4 ) =TT (1,N)-I(N,4):B(1,4) =TT (1,N)-I(N,4)+T(1,1)+ 
0(1,1)-T(N,4)-0(N,4)  ELSE 

A ( 1 , 4 ) =TT ( 1 , N ) -I(N,4)+T(1, 1) +0(1,1) -T(N,4) -0(N,4) :B 
(1,4) =TT (1,N)-I(N,4) 

1060  C ( 1 , 4 ) =1 (N, 4 ) +0 (N , 4 ) +T (N , 4 ) -TT ( 1 , N) 

1070  IF  T ( 1 , 1) +1 (1,1) >T (N,4)+I(N,4)  THEN 

A ( 2 , 4 ) =C-TT (N , 1) : B ( 2 , 4 ) =C-TT (N, 1) +T (1 , 1) +1 (1, 1) -T (N 
, 4 ) -I (N, 4 ) ELSE  A ( 2 , 4 ) = 

C-TT ( N , 1 ) +T ( 1 , 1 ) +1 ( 1 , 1 ) -T (N, 4) -I (N, 4) : B ( 2 , 4 ) =C-TT  (N 

/l) 

1080  C ( 2 , 4 ) =1  (N,  4 ) +T  (N , 4 ) +TT  (N , 1)  -C 

1090  REM  LELA/LALA  ROUTINE 

1100  IF  0(1,3) +T (1,3) >0 ( N , 1 ) +T ( N , 1 ) THEN 

A ( 1 , 9 ) =TT ( 1 , N) : B ( 1 , 9 ) =TT ( 1 , N) +0 ( 1 , 3 ) +T ( 1 , 3 ) -O (N , 1) - 
T (N, 1)  ELSE 

A ( 1 , 9 ) =TT ( 1 , N) +0 (1,3) +T (1,3)-0(N,1) -T (N, 1) : B (1, 9) =T 
T ( 1 , N) 

1110  C ( 1 , 9 ) =0 ( N , 1 ) +T ( N , 1 ) -TT ( 1 , N) 

1120  IF  T ( 1 , 3 ) +1 (1,3) >T (N,1)+I(N,1)  THEN 

A ( 2 , 9 ) =C-TT ( N , 1 ) +0 (1,3) : B ( 2 , 9 ) =C-TT ( N , 1 ) +0 (1,3) +T ( 1 
, 3) +1 (1, 3) -T(N, 1) -I (N, 1)  ELSE 

A ( 2 , 9 ) =C-TT (N , 1 ) +0 (1,3) +T (1,3)+I(1,3) -T (N, 1) -I (N, 1) 

: B (2 , 9 ) =C-TT ( N , 1 ) +0 ( 1 , 3 ) 

1130  C ( 2 , 9 ) =T (N, 1) +1 (N , 1 ) -O (1,3) +TT (N, 1) -C 

1140  REM  LELA/LELE  ROUTINE 

1150  IF  0(1,3) +T (1,3) >0 ( N , 2 ) +T ( N , 2 ) THEN 

A ( 1 , 1 0 ) =TT ( 1 , N ) -I (N, 2 ) : B ( 1 , 1 0 ) =TT ( 1 , N ) -I (N , 2 ) +0 ( 1 , 3 
) +T (1,3) -0 (N , 2 ) -T (N , 2 ) ELSE 
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A ( 1 , 10) =TT (1,N)-I(N,2) +0 (1,3) +T (1,3) -0 (N,  2 ) -T (N,  2 ) : 
B ( 1 , 1 0 ) =TT ( 1 , N ) -I ( N , 2 ) 

1160  C(l, 10)=I(N,2)+O(N,2)+T(N,2) -TT(1,N) 

1170  IF  T ( 1 , 3 ) +1 ( 1 , 3 ) > T(N, 2) +1 (N, 2)  THEN 

A ( 2 , 10) =C-TT (N , 1 ) +0 (1,3) -O (N , 2 ) : B ( 2 , 10 ) =C-TT (N , l)+0 

(1.3)  -0 (N,  2 ) +T (1,3)+I(1,3) -T (N,2)-I(N,2)  ELSE 

A ( 2 , 10 ) =C-TT (N , 1 ) +0 (1,3) -O (N , 2 ) +T (1,3)+I(1,3) -T (N , 2 
) -I (N , 2 ) : B (2 , 10) =C-TT (N , 1) +0 ( 1 , 3 ) -O (N, 2 ) 

1180  C ( 2 , 10) =0 (N , 2 ) +T (N,2)+I(N,2) -O (1,3) +TT (N , 1) -C 

1190  REM  LELA/LALE  ROUTINE 

1200  IF  0(1,3) +T ( 1 , 3 ) >0 ( N , 4 ) +T ( N , 4 ) THEN 

A ( 1 , 12 ) =TT ( 1 , N) -I (N , 4 ) : B ( 1 , 12 ) =TT ( 1 , N) -I (N , 4 ) +0  ( 1 , 3 
) +T (1,3) -T (N,4)-0(N,4)  ELSE 

A ( 1 , 12 ) =TT  ( 1 , N)  -I  (N , 4 ) +0  (1,3)  +T  (1,3)  -0  (N , 4 ) -T  (N , 4 ) : 
B ( 1, 12) =TT ( 1 ,N) -I (N, 4 ) 

1210  C(l, 12)=0(N, 4)+T(N, 4) +1 (N , 4 ) -TT ( 1 , N) 

1220  IF  T ( 1 , 3 ) +1 (1,3) >T (N,4)+I(N,4)  THEN 

A ( 2 , 12 ) =C-TT ( N , 1 ) +0 (1,3) : B ( 2 , 1 2 ) =C-TT (N, 1) +0 ( 1 , 3 ) +T 

(1.3)  +1(1, 3) -T (N,4)-I(N,4)  ELSE 

A ( 2 , 12 ) =C-TT (N , 1) +0 (1,3) +T (1,3) +1 (1,3) -T (N,4)-I(N,4 
) : B ( 2 , 12 ) =C-TT (N , 1 ) +0 (1,3) 

1230  C ( 2 , 12 ) =T  (N,  4 ) +1  (N , 4 ) -0  (1,3)  +TT  (N,  1)  -C 

1240  REM  LELA/LELA  ROUTINE 

1250  IF  (0(1,3)  +T (1,3)  > 0 (N, 3 ) +T (N, 3 ) ) THEN 
A ( 1 , 11)  = 

TT ( 1 , N) : B ( 1 , 11) =TT ( 1 , N) +0 ( 1 , 3 ) +T ( 1 , 3 ) -0 (N , 3 ) -T (N , 3 ) 
ELSE 

A ( 1 , 1 1 ) =TT ( 1 , N ) +0 (1,3) +T (1,3) -O ( N , 3 ) -T (N,3)  :B(1,11) 
=TT(1,N) 

1260  C(l, 11)=0(N, 3) +T(N, 3) -TT(1,N) 

1270  IF  1(1,3) +T (1,3)  > I (N , 3 ) +T (N, 3 ) THEN  A(2,ll)  = 

0(1,3) -0 ( N , 3 ) +C-TT ( N , 1 ) : B (2 , 11)  =T  (1, 3 ) +1  (1 , 3 ) -I  (N,  3 
) -T (N , 3 ) +0 (1,3) -O (N , 3 ) +C-TT (N , 1)  ELSE 
A ( 2 , 11) =T (1,3)+I(1,3) -T (N,3)-I(N,3) +0 (1,3) -0 (N , 3 ) +C 
-TT  (N,  1)  : B ( 2 , 1 1 ) =0  ( 1 , 3 ) -O  ( N , 3 ) +C-TT  ( N , 1 ) 

1280  C ( 2 , 11) =0 (N,3)+I(N,3) +T (N, 3 ) -O (1,3) +TT (N , 1 ) -C 

1290  REM  LALE/LALA  ROUTINE 

1300  IF  0(1,4) +T ( 1 , 4 ) >0 ( N , 1 ) +T ( N , 1 ) THEN 

A ( 1 , 13 ) =TT ( 1 , N) +1 ( 1 , 4 ) : B ( 1 , 13 ) =A ( 1 , 13 ) +0 ( 1 , 4 ) +T  ( 1 , 4 
) -0 (N, 1) -T (N, 1)  ELSE 

B ( 1 , 1 3 ) =TT (1,N)+I(1,4) : A ( 1 , 13 ) =B (1,13) +0 (1,4) +T (1,4 
) -0  (N,  1)  -T  (N,  1) 

1310  C(l, 13)=T(N, 1)+0(N, 1) -1(1,4) -TT ( 1 , N) 

1320  IF  1(1,4) +T (1,4)>I(N,1) +T ( N , 1 ) THEN 

A (2 , 13 ) =C-TT (N, 1)  :B(2, 13)=A(2, 13) +1 (1, 4) +T(1, 4) -I (N 
, 1) -T (N, 1)  ELSE 

B (2 , 13 ) =C-TT (N, 1) : A ( 2 , 13 ) =B (2,13)+I(1,4) +T ( 1 , 4 ) -I (N 
, 1) ~T(N, 1) 

1330  C(2,13)=T(N,1)+I(N, 1)+TT(N, 1) -C 

1340  REM  LALE/LELE  ROUTINE 

1350  IF  0(1,4) +T (1,4) >0 ( N , 2 ) +T ( N , 2 ) THEN 

A ( 1 , 14 ) =TT ( 1 , N) +1 (1,4) -L (N,2)-I(N,2)  : B (1 , 14) =A ( 1 , 14 
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) +0 ( 1 , 4 ) +T ( 1 , 4 ) -0 (N , 2 ) -T ( N , 2 ) ELSE 

B(l, 14)=TT(1,N)+I (1,4) -L(N,2) -I (N, 2 ) : A ( 1 , 14 ) =B ( 1 , 14 
) +0 (1,4) +T (1,4)-0(N,2)-T(N,2) 

1360  C ( 1 , 14 ) =L (N , 2 ) +1 (N , 2 ) +0 (N , 2 ) +T (N,2)-I(l,4) -T (1,4) -T 

T (1 , N) 

1370  IF  I(1,4)+T(1,4)>I(N,2)+T(N,2)  THEN 

A ( 2 , 14 ) =C-TT (N , 1) -L (N, 2 ) -0  (N , 2 ) : B (2 , 14) =A(2 , 14) +1 (1 
,4)+T(l,4)-I(N,2)-T(N,2)  ELSE 

B ( 2 , 14  ) =C-TT  (N , 1)  -L  (N , 2 ) -O  (N , 2 ) : A ( 2 , 14  ) =B  ( 2 , 14 ) +1  ( 1 
, 4 ) +T  (1,4)  -I  (N,  2 ) -T  (N,  2 ) 

1380  C ( 2 , 14 ) =L  (N,  2 ) +0  (N , 2 ) +1 (N , 2 ) +T  (N , 2 ) +TT  (N,  1 ) -C 

1390  REM  LALE/LELA  ROUTINE 

1400  IF  0(1,4) +T ( 1 , 4 ) >0 ( N , 3 ) +T ( N , 3 ) THEN 

A ( 1 , 15 ) =TT ( 1 , N) -I ( 1 , 4 ) : B ( 1 , 15 ) =A ( 1 , 15) +0 ( 1 , 4 ) +T ( 1 , 4 
) -0 (N , 3 ) -T  (N , 3 ) ELSE 

B ( 1 , 15) =TT ( 1 , N ) -1(1,4)  :A(1, 15) =B(1, 15) +0(1,4) +T (1,4 
) -O  (N , 3 ) -T  (N , 3 ) 

1410  C ( 1 , 15 ) =0 (N , 3 ) +T (N,3)-I(l,4) -TT ( 1 , N) 

1420  IF  1(1,4) +T (1,4)>I(N,3) +T ( N , 3 ) THEN 

A ( 2 , 15 ) =C-TT (N , 1) -O (N , 3 ) : B ( 2 , 15 ) =A ( 2 , 15 ) +1 ( 1 , 4 ) +T ( 1 
,4)-I(N,3) -T (N , 3 ) ELSE 

B ( 2 , 1 5 ) =C-TT ( N , 1) -0(N,3) :A(2, 15)=B(2, 15 ) +1 ( 1 , 4 ) +T ( 1 
, 4 ) -I (N , 3 ) -T (N , 3 ) 

1430  C ( 2 , 15 ) =0  (N,  3 ) +1  (N , 3 ) +T  (N , 3 ) +TT  (N,  1)  -C 

1440  REM  LALE/LALE  ROUTINE 

1450  IF  0(1,4) +T (1,4) >0 ( N , 4 ) +T ( N , 4 ) THEN 

A ( 1 , 16 ) =TT ( 1 , N) +1 (1,4)-I(N,4) : B ( 1 , 16) =A ( 1, 16) +0 (1, 4 
)+T(l,4)-0(N,4)-T(N,4)  ELSE 

B ( 1 , 16 ) =TT ( 1 , N) +1 (1,4)-I(N,4) : A ( 1 , 16 ) =B ( 1 , 16 ) +0 ( 1 , 4 
) +T  (1,4) -O (N, 4 ) -T (N, 4 ) 

1460  C ( 1 , 16 ) =1 (N , 4 ) +0 (N , 4 ) +T (N,4)-I(l,4) -TT ( 1 , N) 

1470  IF  1(1,4) +T (1,4)>I(N,4) +T ( N , 4 ) THEN 

A ( 2 , 1 6 ) =C-TT ( N , 1 ) : B (2 , 16 ) =C-TT (N,l)+I(l,4)+T(l,4)-I 
(N , 4 ) -T (N , 4 ) ELSE 

B ( 2 , 16) =C-TT ( N , 1 ) : A ( 2 , 16 ) =C-TT (N,l)+I(l,4) +T ( 1 , 4 ) -I 
(N, 4) -T (N, 4 ) 

1480  C(2, 16) =1 ( N , 4 ) +T ( N , 4 ) +TT ( N , 1) -C 

1540  LPRINT  "MAXBAND  AND  OFFSET  RESULTS" 

1550  LPRINT  TAB (28) /"OFFSET  RANGES  FROM" 

1560  LPRINT  TAB (5) ; "MAX  SUM  OF  BANDS" /TAB (25) ; "OFFSET  LB 
"/TAB (37) /"TO" /TAB (40) "OFFSET  UB" 

1570  FOR  H=1  TO  16 

1580  GOSUB  3830 

1590  LPRINT 

TAB ( 5 ) /MAXB(H) /TAB(25) /OFFB(H) /TAB(40) /OFFE(H) 

1600  NEXT  H 

1660  GOSUB  1730 

1670  FOR  M=2  TO  N-l 

1680  GOSUB  3300 

1690  GOSUB  2510 

1700  GOSUB  3630 

1705  GOSUB  3530 


1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 
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NEXT  M 
END 

OPTBAND=-l 
FOR  E = 1 TO  16 

IF  OPTBAND>  MAXB(E)  THEN  GOTO  1790 
OPTBAND=MAXB ( E ) 

OPTLB  = OFFB(E) 

OPTUP  = OFFE(E) 

NEXT  E 

FOR  E = 1 TO  16 

IF  OPTBAND=MAXB (E)  THEN  GOSUB  1840 
NEXT  E 
RETURN 

LPRINT  "OFFSET  RANGE  FROM 

";OFFB(E) ; "TO" ;OFFE (E) /"YIELDS  MAX  SUM  OF 
BANDS" ;MAXB(E) 

IF  E <=4  THEN  GOSUB  1980: GOTO  1890 
IF  E <=8  THEN  GOSUB  2040: GOTO  1890 
IF  E <=12  THEN  GOSUB  2100:GOTO  1890 
GOSUB  2160 

IF  E MOD  4=1  THEN  GOSUB  2220: GOTO  1930 

IF  E MOD  4=2  THEN  GOSUB  2290: GOTO  1930 

IF  E MOD  4=3  THEN  GOSUB  2360:GOTO  1930 

IF  E MOD  4=0  THEN  GOSUB  2430 

IF  SN0RMAX>S10RMAX  THEN  OCURMAX=SNORMAX  ELSE 
OCURMAX=S 10RMAX 

IF  SNIRMAX>S1IRMAX  THEN  ICURMAX=SNIRMAX  ELSE 
I CURMAX= S 1 I RMAX 

IF  SN0RMIN<S10RMIN  THEN  OCURMIN=SNORMIN  ELSE 
0CURMIN=S10RMIN 

IF  SNIRMIN<S1IRMIN  THEN  ICURMIN=SNIRMIN  ELSE 

ICURMIN=S1IRMIN 

RETURN 

LPRINT  " LAG-LAG  AT  SIGNAL  1"; 

OFFSET (1)=0 

OPTPAT ( 1 ) = l:J=l: GOSUB  3240 

S10RMAX  =OFFSET ( 1 ) : S10RMIN=0FFSET  ( 1) +T  (1, 1) +0  ( 1, 1) 
S1IRMAX=0FFSET ( 1) : S1IRMIN=0FFSET ( 1) +T ( 1 , 1) +1 ( 1 , 1) 
RETURN 

LPRINT  " LEAD-LEAD  AT  SIGNAL  1"; 

OFFSET (1)=0 

OPTPAT ( 1)  = 2:J=l:GOSUB  3240 
S10RMAX 

=OFFSET ( 1 ) +L ( 1 , 2 ) +1 ( 1 , 2 ) : S10RMIN=0FFSET ( 1 ) +L ( 1 , 2 ) +1 

(1.2)  +T (1,2) +0 (1,2) 

S1IRMAX=0FFSET ( 1) +L ( 1 , 2 ) +0 ( 1 , 2 ) : S1IRMIN=0FFSET ( 1) +L 

(1.2)  +0 (1,2) +T (1,2) +1(1, 2) 

RETURN 

LPRINT  " LEAD-LAG  AT  SIGNAL  1"; 

OFFSET (1)=0 

OPTPAT (1)  = 3 : J=1 : GOSUB  3240 

S10RMAX  =OFFSET ( 1 ) : S10RMIN=0FFSET ( 1) +T ( 1 , 3 ) +0 ( 1 , 3 ) 
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S1IRMAX=0FFSET ( 1) +0 ( 1 , 3 ) : S1IRMIN=0FFSET ( 1) +0 ( 1 , 3 ) +T 

(1.3) +I(1,3) 

RETURN 

LPRINT  " LAG-LEAD  AT  SIGNAL  1"; 

OFFSET ( 1) =0 

OPTPAT ( 1)  = 4 : J=1 : GOSUB  3240 
S10RMAX 

=OFFSET ( 1 ) +1 ( 1 , 4 ) : S10RMIN=0FFSET ( 1 ) +1 ( 1 , 4 ) +T ( 1 , 4 ) +0 

(1.4) 

S1IRMAX=0FFSET ( 1 ) : S1IRMIN=0FFSET ( 1 ) +T ( 1 , 4 ) +1 ( 1 , 4 ) 
RETURN 

LPRINT  " LAG-LAG  AT  SIGNAL" ;N 
LPRINT 

OFFSET (N)=OPTLB 

OPTPAT (N)  = l:J=N: GOSUB  3240 

SNORMAX 

=OFFSET (N) -TT ( 1 , N) : SNORMIN=OFFSET (N)  +T  (N,  1)  +0  (N,  1)  - 
TT ( 1 , N ) 

SNIRMAX=OFFSET (N) +TT (N, 1) -C : SNIRMIN=OFFSET (N) +T (N, 1 

) +1 (N, 1) +TT (N, 1) -C 

RETURN 

LPRINT  " LEAD-LEAD  AT  SIGNAL" ;N 
LPRINT 

OFFSET (N)=OPTLB 

OPTPAT (N)  = 2 : J=N : GOSUB  3240 

SNORMAX=OFFSET(N) +L(N,2) +1 (N , 2 ) -TT ( 1 , N) : 

SNORMIN=OFFSET (N) +L (N , 2 ) +1 (N , 2 ) +T (N , 2 ) +0 (N , 2 ) -TT ( 1 , 

N) 

SNIRMAX=OFFSET (N) +L (N , 2 ) +0 (N , 2 ) +TT (N , 1 ) -C : 
SNIRMIN=OFFSET (N ) +L ( N , 2 ) +0 ( N , 2 ) +T ( N , 2 ) +1 (N , 2 ) +TT (N , 
1)-C 
RETURN 

LPRINT  " LEAD-LAG  AT  SIGNAL" ;N 
LPRINT 

OFFSET (N)=OPTLB 

OPTPAT (N)  = 3:J=N:GOSUB  3240 

SNORMAX 

=OFFSET (N) -TT ( 1 , N) : SNORMIN=OFFSET (N) +T (N , 3 ) +0 (N, 3 ) - 
TT  ( 1 , N) 

SNIRMAX=OFFSET (N)  +0 (N , 3 ) +TT (N , 1 ) -C : SNIRMIN=OFFSET (N 

) +0(N,3)+T(N, 3) +1 (N, 3) +TT(N, 1) -C 

RETURN 

LPRINT  " LAG-LEAD  AT  SIGNAL" ;N 
LPRINT 

OFFSET (N)=OPTLB 

OPTPAT (N)  = 4 : J=N : GOSUB  3240 

SNORMAX=OFFSET (N) +1 (N, 4 ) -TT ( 1 , N) : SNORMIN=OFFSET (N) + 
I ( N , 4 ) +T ( N , 4 ) +0 ( N , 4 ) -TT ( 1 , N ) 

SNIRMAX=OFFSET (N) +TT (N , 1 ) -C : SNIRMIN=OFFSET (N) +T (N , 4 

) +1 (N , 4 ) +TT (N , 1 ) -C 

RETURN 

END 
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REM  CALCULATE  MAXIMUM  GREENBAND  AND  OFFSETS 
R=M 

FOR  T = 1 TO  4 

FL  (R, T) =0 : PZ (R, T) =0 : REM  CHECK  FOR  0 BANDWIDTH 
IF  X ( 1 , R, T) <0  THEN 

Z ( 1 , R, T) =X (1,R,T)+Z(1,R,T) :X(1,R,T)=0 
IF  X ( 2 , R,  T) <0  THEN 

Z(2,R,T)=X(2,R,T)+Z(2,R,T) :X(2,R,T)=0 
IF  X(1,R,T)+Z(1,R,T)<0  THEN 
X(1,R,T)=0:Z(1,R,T)=0:Y(1,R,T) =C 
IF  Y (2  , R,  T) <0  THEN 

Z(2,R,T)=Y(2,R,T)+Z(2,R,T) :Y(2,R,T)=0 
IF  X(2,R,T)+Z(2,R,T)<0  THEN 
X (2 ,R, T) =0 : Z (2 , R,T) =0 : Y (2 , R,T) =C 
IF  Y ( 1 , R, T) <0  THEN 

Z ( 1 , R, T) =Y (1,R,T)+Z(1,R,T) : Y ( 1 , R, T) =0 

IF  Z(1,R,T)<0  THEN  W ( 1 , R, T) =-Z ( 1 , R, T)  ELSE 
W ( 1 , R , T) =0 

IF  X(1,R,T)+Y(1,R,T)+Z(1,R,T)  < C THEN 
Q ( 1 , R, T) =X ( 1 , R, T) +Y (1,R,T)+Z(1,R,T)  ELSE 
Q(1,R,T)=C 

IF  Z (2 , R, T) <0  THEN  W ( 2 , R, T) =-Z ( 2 , R, T)  ELSE 
W (2 , R, T) =0 

IF  X ( 2 , R, T) +Y (2,R,T)+Z(2,R,T)  < C THEN 
Q ( 2 , R, T) =X ( 2 , R, T) +Y (2,R,T)+Z(2,R,T)  ELSE 
Q(2,R,T)=C 

IF  W ( 2 , R,  T) <=  W ( 1 , R, T)  THEN 
PZ (R, T) =PZ (R,T)+1000  ELSE  IF  W ( 2 , R, T) <X ( 1 , R, T) 

THEN  PZ (R, T) =PZ (R,T)+2000  ELSE  IF 

W (2 , R, T) <Y ( 1 , R, T)  THEN  PZ (R, T) =PZ (R, T) +3000  ELSE 
IF  W (2 , R, T) <Q ( 1 , R, T)  THEN  PZ (R, T) =PZ (R, T) +4000 
ELSE  PZ (R, T) =PZ (R, T) +5000 

IF  X ( 2 , R, T) < W ( 1 , R, T)  THEN  PZ (R, T) =PZ (R, T) +100  ELSE 
IF  X(2,R,T)<X(1,R,T)  THEN  PZ (R, T) =PZ (R, T) +200  ELSE 
IF  X(2,R,T)<Y(1,R,T)  THEN  PZ (R, T) =PZ (R, T) +3 00  ELSE 
IF  X(2,R,T)<Q(1,R,T)  THEN  PZ (R, T) =PZ (R, T) +400  ELSE 
PZ (R, T) =PZ (R, T) +500 

IF  Y ( 2 , R, T) < W ( 1 , R, T)  THEN  PZ (R, T) =PZ (R, T) +10 
ELSE  IF  Y(2,R,T)<X(1,R,T)  THEN  PZ (R, T) =PZ (R, T) +20 
ELSE  IF  Y(2,R,T)<Y(1,R,T)  THEN  PZ (R, T) =PZ (R, T) +30 
ELSE  IF  Y (2 , R, T) <Q ( 1 , R, T)  THEN  PZ (R, T) =PZ (R, T) +4 0 
ELSE  PZ  (R,  T)  =PZ  (R,  T)  +50 

IF  Q ( 2 , R,  T)  < W ( 1 , R,  T)  THEN  PZ  (R,  T)  =PZ  (R,  T) +1 
ELSE  IF  Q (2  , R,  T)  <X  ( 1 , R,  T)  THEN  PZ  (R,  T)  =PZ  (R,  T)  +2 
ELSE  IF  Q(2,R,T)<Y(1,R,T)  THEN  PZ (R, T) =PZ (R, T) +3 
ELSE  IF  Q(2,R,T)<Q(1,R,T)  THEN  PZ (R, T) =PZ (R, T) +4 
ELSE  PZ (R, T) =PZ (R, T) +5 
IF  PZ (R, T) =1111  THEN  IF 

X(1,R,T)+Z(1,R,T) >X (2,R,T)+Z(2,R,T)  THEN 

GRBD  (R,  T)  =X  ( 1 , R,  T)  +Z  ( 1 , R,  T)  : VARB  (R,  T)  =X  ( 1 , R,  T)  : VARE 

(R,  T)  =Y  ( 1 , R,  T)  ELSE 
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GRBD (R , T) =X (2,R,T)+Z(2,R,T) : VARB (R, T) =X ( 2 , R, T) : VARE 
(R,T)=Y(2,R,T) 

IF  PZ (R,  T) =1112  THEN  GRBD (R, T) = 

X(1,R,T)+Z(1,R,T) : VARB ( R , T ) =X ( 1 , R , T ) : 

VARE ( R , T ) = Y ( 1 , R , T ) 

IF  PZ (R, T) =1113  THEN  GRBD (R, T) = 

X(2,R,T)+Z(2,R,T) : VARB (R, T) =X ( 2 , R, T) : 

VARE (R,T)=Y ( 2 , R , T ) 

IF  PZ (R, T) =1114  THEN  GRBD (R, T) = 

X (2  , R,  T)  +Z  ( 2 , R, T)  : VARB (R, T) =X ( 2 , R,  T)  : 

VARE ( R , T ) = Y ( 2 , R , T ) 

IF  PZ (R,  T) =1115  THEN  GRBD (R, T) = 

X(2,R,T)+Z(2,R,T) : VARB (R, T) =X ( 2 , R, T) : 

VARE ( R , T ) = Y ( 2 , R , T ) 

IF  PZ (R, T) =1122  THEN  GRBD (R, T) = 

X(1,R,T)+Z(1,R,T) : VARB ( R , T ) =X ( 1 , R , T ) : VARE (R, T) =Y ( 1 , 
R,T) 

IF  PZ (R, T) =1123  THEN  GRBD (R, T) = 

X (2  , R,  T) +Y (2,R,T)+Z(2,R,T)+Z(1,R,T) : VARB (R, T) =Y ( 2 , R 
, T) : VARE (R,T)=X(1,R,T) 

IF  PZ (R,  T) =1124  THEN  GRBD (R, T) = 

X(2,R,T)+Y(2,R,T)+Z(2,R,T)+Z(1,R,T) : VARB (R, T) =Y ( 2 , R 
, T) : VARE ( R , T ) =X ( 1 , R , T ) 

IF  PZ (R, T) =1125  THEN  GRBD (R, T) = 

X(2,R,T)+Y(2,R,T)+Z(2,RfT)+Z(l,R,T) : VARB (R, T) =Y ( 2 , R 
, T) : VARE ( R , T ) =X ( 1 , R , T ) 

IF  PZ (R,T) =1133  THEN  GRBD (R; T) = 

X(2,R,T)+X(1,R,T)+Z(2,R,T)+Z(1,R,T) : VARB (R, T) =X ( 1 , R 
, T) : VARE ( R , T ) = Y ( 2 , R , T ) 

IF  PZ (R, T) =1134  THEN  GRBD (R, T) = 

X(2,R,T) +X(1,R,T)+Z (2,R,T)+Z (1,R,T) : VARB (R, T) =X ( 1 , R 
, T) : VARE ( R , T ) = Y ( 2 , R , T ) 

IF  PZ (R,T) =1135  THEN  GRBD (R, T) = 

X(2fR,T)+X(l,RfT)+Z(2,R,T)+Z(l,R,T) : VARB (R, T) =X ( 1 , R 
, T) : VARE (R,T)=Y(1,R,T) 

IF  PZ (R, T) =1144  THEN  GRBD (R, T) = 

X (2 , R, T) +X (1/R,T)+Z(2,R,T)+Z(1,R,T) : VARB (R,T) =X(1,R 
, T) : VARE (R,T)=Y(1,R,T) 

IF  PZ (R, T) =1145  THEN  GRBD (R, T) = 

X ( 2 , R, T) +X (1,R,T)+Z(2,R,T)+Z(1,R,T) : VARB (R, T) =X ( 1 , R 
, T) : VARE ( R , T ) =Y ( 1 , R , T ) 

IF  PZ (R, T) =1155  THEN  GRBD (R, T) = 

X (2 , R, T) +X (1,R,T)+Z(2,R,T)+Z(1,R,T) : VARB (R, T) =X ( 1 , R 
, T) : VARE (R, T) =Y ( 1 , R, T) 

IF  PZ (R, T) =1222  THEN 

GRBD (R, T) =X ( 1 , R, T) +Z(1,R,T) : VARB (R, T) =X ( 1 , R, T) : VARE 
(R/  T)  =Y  ( 1 , R,  T) 

IF  PZ (R, T) =1223  THEN 

GRBD ( R , T ) =X ( 2 , R , T ) +Y ( 2 , R, T) +Z ( 2 , R, T) +Z ( 1 , R, T) :VARB( 
R, T) =Y (2 , R, T) : VARE ( R , T ) =X ( 1 , R , T ) 
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IF  PZ (R,  T) =1224  THEN 

GRBD (R, T) =X (2 , R, T) +Y (2,R,T)+Z(2,R,T)+Z(1,R,T) :VARB( 
R, T) =Y (2 , R,  T) : VARE ( R , T ) =X ( 1 , R , T ) 

IF  PZ (R, T) =1225  THEN 

GRBD (R,  T) =X (2 , R, T) +Y (2,R/T)+Z(2,R,T)+Z(1/R,T) :VARB( 
R, T) =Y (2 , R, T) : VARE ( R , T ) =X ( 1 , R , T ) 

IF  PZ (R, T) =1233  THEN 

GRBD (R, T) =X (2,R,T)+Z(2,R,T) +X (1,R,T)+Z(1,R,T) :VARB( 
R,  T)  =X  ( 1 , R,  T)  : VARE  ( R , T ) = Y ( 2 , R , T ) 

IF  PZ (R, T) =1235  THEN 

GRBD (R,  T) =X (2,R,T)+Z(2,R,T) +X (1,R,T)+Z(1,R,T) :VARB( 
R,  T) =X ( 1 , R, T) : VARE ( R , T ) =Y ( 2 , R , T ) 

IF  PZ (R,  T) =1244  THEN 

GRBD (R,  T) =X (2,R,T)+Z(2,R,T) +X ( 1 , R, T) +Z ( 1 , R, T) :VARB( 
R, T) =X ( 1 , R, T) : VARE ( R , T ) = Y ( 1 , R , T ) 

IF  PZ (R, T) =1245  THEN 

GRBD (R,  T) =X (2,R,T)+Z(2,R,T) +X (1,R,T)+Z(1,R,T) :VARB( 
R, T) =X ( 1 , R, T) : VARE ( R , T ) = Y ( 1 , R , T ) 

IF  PZ (R,  T) =1255  THEN 

GRBD ( R , T ) =X ( 2 , R , T ) +Z ( 2 , R , T) +X ( 1 , R, T) +Z ( 1 , R, T) :VARB( 
R, T) =X ( 1 , R, T) : VARE ( R , T ) = Y ( 1 , R , T ) 

IF  PZ (R, T) =1333  THEN 

GRBD (R, T) =X (2,R,T)+Z(2,R,T) +X (1,R,T)+Z(1,R,T) :VARB( 
R,  T)  =X  (2  , R, T)  : VARE  ( R , T ) = Y ( 2 , R , T ) 

IF  PZ (R, T) =1334  THEN 

GRBD (R, T) =X (2,R,T)+Z(2,R,T) +X ( 1 , R, T) +Z ( 1 , R, T) :VARB( 
R, T) =X (2 , R, T) : VARE (R,T)=Y(2,R,T) 

IF  PZ (R,T) =1335  THEN 

GRBD (R, T) =X (2,R,T)+Z(2,R,T) +X ( 1 , R, T) +Z ( 1 , R, T) :VARB( 
R,  T)  =X  ( 2 , R,  T)  : VARE  ( R , T ) = Y ( 2 , R , T ) 

IF  PZ (R, T) =1344  THEN 

GRBD (R, T) =X (2,R,T)+Z(2,R,T) +X (1,R,T)+Z(1/R,T) :VARB( 
R, T) =X ( 2 , R, T) : VARE ( R , T ) =Y ( 1 , R , T ) 

IF  PZ (R, T) =1345  THEN 

GRBD (R, T) =X (2,R,T)+Z(2,R,T) +X (1/R,T)+Z(1,R,T) :VARB( 
R,  T)  =X  ( 2 , R,  T)  : VARE  ( R , T ) = Y ( 1 , R , T ) 

IF  PZ (R, T) =1355  THEN 

GRBD (R, T) =X (2 , R, T) +Z (2 , R, T) +X (1,R,T)+Z(1,R,T) :VARB( 
R, T) =X ( 2 , R, T)  : VARE ( R , T ) = Y ( 1 , R , T ) 

IF  PZ (R, T) =1444  THEN 

GRBD (R, T) =X ( 1 , R, T) +Y (1,R,T)+Z(1,R,T)+Z(2/R,T) :VARB( 
R, T) =Y ( 1 , R, T) : VARE ( R , T ) =X ( 2 , R , T ) 

IF  PZ (R, T) =1445  THEN 

GRBD (R, T) =X ( 1 , R, T) +Y (1,R,T)+Z(1,R,T)+Z(2,R,T) :VARB( 
R, T) =Y ( 1 , R, T) : VARE ( R , T ) =X ( 2 , R , T ) 

IF  PZ (R/ T) =1455  THEN 

GRBD (R, T) =X ( 1 , R, T) +Y (1,R,T)+Z(1,R,T)+Z(2,R,T) :VARB( 
R, T) =Y ( 1 , R, T) : VARE (R,T)=X(2,R,T) 

IF  PZ (R, T) =1555  THEN 

GRBD (R, T) =X (2,R,T)+Z(2,R,T) : VARB (R, T) =X (2 , R, T) : VARE 
(R,T)=Y(2,R,T) 
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IF  PZ (R,  T) >1555  THEN  FL (R, T) =9 : GOSUB  4360  ELSE  GOTO 
3050 

IF  FL (R, T) =9  THEN  GOTO  2510 

NEXT  T 

REM  RETURN 

RETURN 

HV1H=V ( 1 , H) 

HA1H=A ( 1 , H) 

HB1H=B ( 1 , H) 

HU1H=U ( 1 , H) 

HC1H=C ( 1 , H) 

V(1,H)=V(2,H) 

A ( 1 / H)  =A  ( 2 , H) 

B ( 1 / H) =B (2 , H) 

U ( 1 / H) =U ( 2 , H) 

C(1,H)=C(2,H) 

V ( 2 , H) =HV1H 
A ( 2 , H) =HA1H 
B ( 2 , H) =HB1H 
U (2 , H) =HU1H 
C ( 2 , H) =HC1H 
RETURN 

REM  COMPUTE  BEG  AND  END  OF  BANDS 
IF  OPTPAT ( J) =1  THEN 

BOB ( J) =OFFSET ( J) : EOB ( J) =OFFSET ( J) +T(J, 1) +0(J, 1) :BIB 
( J) =OFFSET ( J) : EIB ( J) =OFFSET ( J) +T ( J; 1) +1 ( J, 1) 

IF  OPTPAT (J) =2  THEN 

BOB ( J) =OFFSET ( J) +L ( J , 2 ) +1 ( J , 2 ) : EOB ( J) =OFFSET ( J) +L ( J 
, 2 ) +1 ( J , 2 ) +0 ( J , 2 ) +T ( J , 2 ) : BIB ( J) =OFFSET ( J) +L ( J, 2) +0 ( 
J , 2 ) : EIB ( J) =OFFSET ( J) +L ( J, 2 ) +1 ( J, 2 ) +0 ( J, 2 ) +T ( J, 2 ) 

IF  OPTPAT (J) =3  THEN 

BOB ( J) =OFFSET ( J) : EOB ( J) =OFFSET ( J) +0 ( J, 3 ) +T ( J, 3 ) :BIB 
( J) =OFFSET ( J) +0(J, 3) : EIB ( J) =OFFSET ( J) +1 ( J, 3 ) +0 ( J, 3 ) 
+T(J,3) 

IF  OPTPAT (J) =4  THEN 

BOB ( J) =OFFSET ( J) +1 (J, 4) : EOB ( J) =OFFSET ( J) +1 ( J, 4 ) +0 ( J 
, 4 ) +T ( J , 4 ) : BIB ( J) =OFFSET ( J) : EIB ( J) =OFFSET ( J) +1 ( J , 4 ) 
+T(J,4) 

RETURN 
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3300  REM  LAG-LAG  ON  MIDDLE  SIGNAL 

3310  IF  T (M, 1) +0 (M, 1) < OCURMIN-OCURMAX  THEN 

X ( 1 , M, 1) =TT ( 1 , M) +OCURMAX: Y(1,M, 1)=TT(1,M) +OCURMIN-T 
(M, 1) -O (M, 1)  ELSE 

X ( 1 ,M, 1) =TT ( 1 ,M) +OCURMIN-T (M, 1)-0(M,1) : Y ( 1 , M , 1 ) =TT ( 

1 ,  M) +OCURMAX 

3320  Z (1,M, 1)=T(M, 1) +0(M, 1) -TT(1,M) -OCURMAX 
3330  IF  I (M, 1) +T (M, 1) <ICURMIN-ICURMAX  THEN 

X ( 2 , M,  1) =C-TT (N , 1 ) +TT (N , M) +ICURMAX : Y (2 , M, 1)=C-TT(N, 

1)  +TT (N,M) +ICURMIN-I (M, 1) -T(M, 1)  ELSE 

X (2 ,M, 1) =C-TT (N, 1) +TT (N,M) +ICURMIN-I (M, 1) -T (M, 1) : Y ( 

2 , M, 1) =C-TT (N, 1) +TT (N, M) +ICURMAX 

3340  Z ( 2 , M, 1) =T (M,1)+I(M,1) -TT (N, M) +TT (N, 1) -C-ICURMAX 

3350  REM  LEAD-LEAD  ON  MIDDLE  SIGNAL 

3360  IF  T (M, 2 ) +0 (M, 2 ) < OCURMIN-OCURMAX  THEN 

X ( 1 , M, 2 ) =TT ( 1 , M) +OCURMAX-L (M,  2 ) -I ( M , 2 ) : Y ( 1 , M , 2 ) =TT ( 

1 , M) +OCURMIN-L (M,2)-I(M,2) -T (M,  2 ) -O (M, 2 ) ELSE 

X ( 1 , M , 2 ) =TT ( 1 , M) +OCURMIN-L (M,2)-I(M,2) -O (M, 2 ) -T (M, 2 
) : Y ( 1 , M, 2 ) =TT ( 1 , M) +OCURMAX-L (M,  2 ) -I (M, 2 ) 

3 370  Z ( 1 , M, 2 ) =L  (M,  2 ) +1 (M, 2 ) +T  (M,  1 ) +0  (M,  1)  -TT  ( 1 , M)  -OCURMA 
X 

3380  IF  I (M, 2 ) +T (M, 2 ) <ICURMIN-ICURMAX  THEN 

X ( 2 , M , 2 ) =C-TT (N, 1) +TT (N , M) +ICURMAX-L (M, 2 ) -O (M, 2 ) :Y( 

2 ,  M , 2 ) =C-TT (N , 1) +TT (N , M) +ICURMIN-L (M, 2 ) -O (M,2)-I(M, 

2 )  -T (M, 2 ) ELSE 

X ( 2 , M, 2 ) =C-TT (N , 1 ) +TT (N , M) +ICURMIN-L (M , 2 ) -O (M, 2 ) -I ( 
M, 2 ) -T (M,  2 ) : Y ( 2 , M , 2 ) =C-TT ( N , 1) +TT (N,M) +ICURMAX- 
L (M, 2 ) -O (M, 2 ) 

3390  Z ( 2 , M,  2 ) =L  (M,  2 ) +0  (M,  2 ) +1  (M,  2 ) +T  (M,  2 ) -TT  (N, M)  +TT  (N,  1 
) -C-ICURMAX 

3400  REM  LEAD-LAG  ON  MIDDLE  SIGNAL 

3410  IF  T (M, 3 ) +0 (M, 3 ) < OCURMIN-OCURMAX  THEN 

X ( 1 ,M, 3 ) =TT ( 1 , M) +OCURMAX : Y ( 1 , M, 3 ) =TT ( 1 , M) +OCURMIN-T 
(M, 3 ) -O (M, 3 ) ELSE 

X ( 1 , M , 3 ) =TT ( 1 , M) +0CURMIN-0 (M,3)-T(M,3):Y(1,M,3) =TT ( 
1 , M) +OCURMAX 

3420  Z ( 1 , M, 3 ) =0 (M, 3 ) +T (M, 3 ) -TT ( 1 , M) -OCURMAX 
3430  IF  I (M, 3 ) +T (M, 3 ) <ICURMIN-ICURMAX  THEN 

X ( 2 , M , 3 ) =C-TT ( N , 1 ) +TT ( N , M ) + I CURMAX-0 (M,3):Y(2,M,3)= 
C-TT (N, 1) +TT(N,M) +ICURMIN-0 (M, 3) -I (M, 3 ) -T (M, 3 ) ELSE 
X (2  , M,  3 ) =C-TT  (N,  1)  +TT  (N,  M)  +ICURMIN-0  (M,  3 ) -I  (M,  3 ) -T  ( 
M,  3 ) : Y ( 2 , M,  3 ) =C-TT  (N , 1)  +TT  (N,  M)  +ICURMAX-0  (M,  3 ) 

3440  Z (2  , M,  3 ) =0  (M,  3 ) +1  (M,  3 ) +T  (M,  3 ) -TT  (N , M) +TT  (N,  1)  -C-ICU 
RMAX 

3450  REM  LAG-LEAD  ON  MIDDLE  SIGNAL 

3460  IF  T (M, 4 ) +0 (M, 4 ) < OCURMIN-OCURMAX  THEN 

X ( 1 , M, 4 ) =TT ( 1 , M) +OCURMAX-I (M,4):Y(1,M,4) =TT ( 1 , M) +OC 
URMIN-T (M, 4 ) -0 (M,4)-I(M,4)  ELSE 

X ( 1 , M, 4 ) =TT ( 1 , M) +0CURMIN-0 (M, 4 ) -T (M, 4 ) -I (M, 4 ) : Y ( 1 , M 
, 4 ) =TT ( 1 , M) +OCURMAX-I (M , 4 ) 

3470  Z ( 1 , M, 4 ) =0 (M, 4 ) +T (M,4)+I(M,4) -TT ( 1 , M) -OCURMAX 
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3480  IF  I (M, 4 ) +T (M, 4 ) <ICURMIN-ICURMAX  THEN 

X ( 2 , M , 4 ) =C-TT (N , 1 ) +TT (N , M) +ICURMAX : Y ( 2 , M, 4 ) =C-TT (N , 
1 ) +TT (N , M) +ICURMIN-I (M, 4 ) -T (M, 4 ) ELSE 
X ( 2 , M , 4 ) =C-TT (N , 1 ) +TT (N , M) +ICURMIN-I (M , 4 ) -T (M , 4 ) : Y ( 
2 , M, 4 ) =C-TT  (N , 1)  +TT  (N , M)  +ICURMAX 
3490  Z (2 , M; 4 ) =1 (M,  4 ) +T (M,  4 ) -TT (N,  M) +TT (N, 1) -C-ICURMAX 

3500  FOR  J=1  TO  4 

3501  IF  X(1,M,J)>C  THEN  X ( 1 , M , J) =X ( 1 , M, J)  MOD 
C:Y(1,M,J)=Y(1,M,J)  MOD  C: Z ( 1M, J) =Z ( 1 , M, J) +C 

3502  IF  X (2 , M, J) >C  THEN  X ( 2 , M, J) =X (2 , M, J)  MOD 

C : Y ( 2 , M , J ) = Y ( 2 , M , J ) MOD  C: Z (2 ,M, J) =Z (2 , M, J) +C 

3503  IF  Y ( 1 , M,  J) >C  THEN  Y(1,M,J)=C 

3504  IF  Y ( 2 , M, J) >C  THEN  Y(2,M,J)=C 

3510  NEXT  J 

3520  RETURN 

3530  REM  RESET  BEG  AND  END  OF  BANDS  MOD  CYCLE  LENGTH 
3550  IF  EOB(M) -TT(1,M)  <0  THEN 

EOB (M) =EOB (M) +C: BOB (M) =BOB (M) +C : GOTO  3550 
3570  IF  EIB(M) -TT  (N , M) +TT  (N , 1) -C  <0  THEN 

EIB (M) =EIB (M) +C : BIB (M) =BIB (M) +C: GOTO  3570 
3580  IF  BOB(M) -TT(1,M)  > OCURMAX  THEN 

OCURMAX=BOB (M) -TT ( 1 , M) 

3590  IF  EOB (M) -TT ( 1 , M)  < OCURMIN  THEN 

OCURMIN=EOB (M) -TT ( 1 , M) 

3600  IF  BIB (M) -TT (N,M) +TT (N, 1) -C  > ICURMAX  THEN 
ICURMAX=BIB(M) -TT (N , M) +TT (N, 1) -C 
3610  IF  EIB(M)-TT(N,M)+TT(N,1)-C  < ICURMIN  THEN 

ICURMIN=EIB (M) -TT (N , M) +TT (N, 1) -C 
3620  RETURN 

3630  REM  DETERMINE  MAX  ON  MID  SIGNALS 
3640  OPTBAND=-l 

3650  FOR  E = 1 TO  4 

3660  IF  OPTBAND>  GRBD (M, E)  THEN  GOTO  3700 

3670  OPTBAND=GRBD(M,E) 

3680  OPTLB  = VARB (M,  E) 

3690  OPTUP  = VARE (M, E) 

3700  NEXT  E 

3710  FOR  E = 1 TO  4 

3720  IF  OPTBAND=GRBD (M, E)  THEN  GOSUB  3750 

3730  NEXT  E 

3740  RETURN 

3750  LPRINT  "OFFSET  RANGE  FROM 

";VARB(M,E) ; "TO"; VARE (M,E) /"YIELDS  MAX  SUM  OF 
BANDS" ; GRBD (M,E) 

3760  IF  E =1  THEN  LPRINT  "LAG-LAG  AT 

SIGNAL" ;M: OFFSET (M) =VARB (M , E) : OPTPAT (M) =1 
3770  IF  E =2  THEN  LPRINT  "LEAD-LEAD  AT 

SIGNAL" ;M: OFFSET (M) =VARB (M, E) : OPTPAT (M) =2 
3780  IF  E =3  THEN  LPRINT  "LEAD-LAG  AT 

SIGNAL" ;M: OFFSET (M) =VARB (M, E) : OPTPAT (M) =3 
3790  IF  E =4  THEN  LPRINT  "LAG-LEAD  AT 

SIGNAL" ;M: OFFSET (M) =VARB (M, E) : OPTPAT (M) =4 


89 


3800  J=M:GOSUB  3240 

3810  REM  GOSUB  3530 

3820  RETURN 

3830  REM  CALCULATE  MAXIMUM  GREENBAND  AND  OFFSETS 
3840  FA (H) =0 : PS(H)=0:REM  CHECK  FOR  0 BANDWIDTH 
3850  IF  A ( 1 , H) <0  THEN  A(1,H)=0 

3860  IF  A ( 2 , H) <0  THEN  A(2,H)=0 

3870  IF  A ( 1 , H) +C ( 1 , H) <0  THEN  A (1, H) =0 : C ( 1 ,H) =0 : B ( 1,H) =C 
3880  IF  B ( 2 , H) <0  THEN  B(2,H)=0 

3890  IF  A ( 2 , H) +C ( 2 , H) <0  THEN  A ( 2 , H) =0 : C ( 2 , H) =0 : B ( 2 , H) =C 
3900  IF  B ( 1 , H) <0  THEN  B(1,H)=0 

3910  IF  C ( 1 , H) <0  THEN  V ( 1 , H) =-C ( 1 , H)  ELSE  V(1,H)=0 

3920  IF  A ( 1 , H) +B ( 1 , H) +C ( 1 , H)  < C THEN 

U ( 1 , H ) =A ( 1 , H ) +B ( 1 , H ) +C ( 1 , H ) ELSE  U(1,H)=C 
3930  IF  C ( 2 , H) <0  THEN  V ( 2 , H) =-C ( 2 , H)  ELSE  V(2,H)=0 
3940  IF  A ( 2 , H) +B (2 , H) +C ( 2 , H)  < C THEN 

U ( 2 , H) =A ( 2 , H) +B ( 2 , H) +C ( 2 , H)  ELSE  U(2,H)=C 
3950  IF  V ( 2 , H) <=  V(1,H)  THEN  PS (H) =PS (H) +1000  ELSE  IF 
V(2,H)<A(1,H)  THEN  PS (H) =PS (H) +2 000  ELSE  IF 
V ( 2 , H) <B ( 1 , H)  THEN  PS (H) =PS (H) +3 000  ELSE  IF 
V ( 2 , H) <U ( 1 , H)  THEN  PS (H) =PS (H) +4 000  ELSE 
PS (H) =PS (H) +5000 

3960  IF  A ( 2 , H) < V(1,H)  THEN  PS (H) =PS (H) +100  ELSE  IF 
A (2 , H) <A ( 1 , H)  THEN  PS (H) =PS (H) +200  ELSE  IF 
A ( 2 , H) <B ( 1 , H)  THEN  PS (H) =PS (H) +300  ELSE  IF 
A ( 2 , H) <U ( 1 , H)  THEN  PS (H) =PS (H) +4 00  ELSE 
PS (H) =PS (H) +500 

3970  IF  B ( 2 , H) < V(1,H)  THEN  PS (H) =PS (H) +10  ELSE  IF 
B ( 2 , H) <A ( 1 , H)  THEN  PS (H) =PS (H) +2 0 ELSE  IF 
B ( 2 , H) <B ( 1 , H)  THEN  PS (H) =PS (H) +3 0 ELSE  IF 
B ( 2 , H) <U ( 1 , H)  THEN  PS (H) =PS (H) +4 0 ELSE 
PS (H) =PS (H) +50 

3980  IF  U(2,H)<  V(1,H)  THEN  PS (H) =PS (H) +1  ELSE  IF 
U ( 2 , H) <A ( 1 , H)  THEN  PS (H) =PS (H) +2  ELSE  IF 
U ( 2 , H) <B ( 1 , H)  THEN  PS (H) =PS (H) +3  ELSE  IF 
U ( 2 , H) <U ( 1 , H)  THEN  PS (H) =PS (H) +4  ELSE  PS (H) =PS (H) +5 
3990  IF  PS (H) =1111  THEN  IF  A ( 1 , H) +C  (1,  H)  >A  (2  , H) +C  (2  , H) 
THEN  MAXB (H) = 

A ( 1 , H) +C ( 1 , H) : OFFB (H) =A ( 1 , H) : OFFE (H) =B ( 1 , H)  ELSE 
MAXB ( H ) =A ( 2 , H ) +C ( 2 , H ) : OFFB (H) =A ( 2 , H) : OFFE (H) =B ( 2 , H) 
4000  IF  PS (H) =1112  THEN  MAXB(H)= 

A ( 1 , H) +C ( 1 , H) : OFFB (H) =A ( 1 , H) : OFFE (H) =B ( 1 , H) 

4010  IF  PS (H) =1113  THEN  MAXB(H)= 

A ( 2 , H) +C ( 2 , H) : OFFB (H) =A ( 2 , H) : OFFE (H) =B ( 2 , H) 

4020  IF  PS (H) =1114  THEN  MAXB(H)= 

A (2 , H) +C ( 2 , H) : OFFB (H) =A ( 2 , H) : OFFE (H) =B ( 2 , H) 

4030  IF  PS (H) =1115  THEN  MAXB(H)= 

A ( 2 , H) +C ( 2 , H) : OFFB (H) =A ( 2 , H) : OFFE (H) =B ( 2 , H) 

4040  IF  PS (H) =1122  THEN  MAXB(H)= 

A(1,H)+C(1,H) : OFFB (H) =A ( 1 , H) : OFFE (H) =B ( 1 , H) 


4050 

4060 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 
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IF  PS (H) =1123  THEN  MAXB(H)= 

A ( 2 , H) +B ( 2 , H) +C ( 2 , H) +C ( 1 , H) : OFFB (H) =B (2 , H) :OFFE(H)= 
A(1,H) 

IF  PS (H) =1124  THEN  MAXB(H)= 

A ( 2 , H) +B ( 2 , H) +C ( 2 , H) +C ( 1 , H) : OFFB (H) =B ( 2 , H) :OFFE(H)= 
A ( 1 / H) 

IF  PS (H) =1125  THEN  MAXB(H)= 

A ( 2 , H) +B ( 2 , H) +C ( 2 , H) +C ( 1 , H) : OFFB (H) =B (2 , H) :OFFE(H)= 
A ( 1 , H) 

IF  PS (H) =1133  THEN  MAXB(H)= 

A ( 2 , H ) +A ( 1 , H ) +C ( 2 , H ) +C ( 1 , H ) : OFFB (H) =A ( 1 , H) :OFFE(H)= 
B (2 , H) 

IF  PS (H) =1134  THEN  MAXB(H)= 

A ( 2 , H) +A ( 1 , H) +C (2 , H) +C ( 1 , H) : OFFB (H) =A ( 1 , H) : OFFE (H) = 
B ( 2 , H) 

IF  PS (H) =1135  THEN  MAXB(H)= 

A ( 2 , H) +A ( 1 , H) +C ( 2 , H) +C ( 1 , H) : OFFB (H) =A ( 1 , H) :OFFE(H)= 
B ( 1 / H) 

IF  PS (H) =1144  THEN  MAXB (H) = 

A ( 2 , H) +A ( 1 , H) +C ( 2 , H ) +C ( 1 , H) : OFFB (H) =A ( 1 , H) :OFFE(H)= 
B(1,H) 

IF  PS (H) =1145  THEN  MAXB(H)= 

A ( 2 , H ) +A ( 1 , H ) +C ( 2 , H ) +C ( 1 , H ) : OFFB (H) =A ( 1 , H) :OFFE(H)= 
B ( 1 / H) 

IF  PS (H) =1155  THEN  MAXB (H) = 

A ( 2 , H ) +A ( 1 , H ) +C ( 2 , H ) +C ( 1 , H ) : OFFB (H) =A ( 1 , H) :OFFE(H)= 
B(1,H) 

IF  PS (H) =1222  THEN 

MAXB ( H ) =A ( 1 , H ) +C ( 1 , H ) : OFFB (H) =A ( 1 , H) : OFFE (H) =B ( 1 , H) 
IF  PS (H) =1223  THEN 

MAXB ( H ) =A ( 2 , H ) +B ( 2 , H ) +C ( 2 , H ) +C ( 1 , H ) : OFFB (H) =B ( 2 , H) : 
OFFE (H) =A ( 1 , H) 

IF  PS (H) =1224  THEN 

MAXB ( H ) =A ( 2 , H ) +B (2,H)+C(2,H) +C ( 1 , H ) : OFFB (H) =B ( 2 , H) : 
OFFE ( H ) =A ( 1 , H ) 

IF  PS (H) =1225  THEN 

MAXB ( H ) =A ( 2 , H ) +B ( 2 , H ) +C ( 2 , H ) +C ( 1 , H ) : OFFB (H) =B ( 2 , H) : 
OFFE (H) =A ( 1 , H) 

IF  PS (H) =1233  THEN 

MAXB ( H ) =A ( 2 , H ) +C ( 2 , H ) +A ( 1 , H ) +C ( 1 , H ) : OFFB (H) =A ( 1 , H) : 
OFFE (H) =B ( 2 , H) 

IF  PS (H) =12  3 5 THEN 

MAXB ( H ) =A ( 2 , H ) +C ( 2 , H ) +A ( 1 , H ) +C ( 1 , H ) : OFFB (H) =A ( 1 , H) : 
OFFE (H) =B ( 2 , H) 

IF  PS (H) =1244  THEN 

MAXB ( H ) =A ( 2 , H ) +C ( 2 , H ) +A (1,H)+C(1,H) : OFFB (H) =A ( 1 , H) : 
OFFE (H) =B ( 1 , H) 

IF  PS (H) =1245  THEN 

MAXB ( H ) =A ( 2 , H ) +C ( 2 , H ) +A ( 1 , H ) +C ( 1 , H ) : OFFB (H) =A ( 1 , H) : 
OFFE (H) =B ( 1 , H) 
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4220  IF  PS (H) =1255  THEN 

MAXB ( H ) =A ( 2 , H ) +C ( 2 , H ) +A ( 1 , H ) +C ( 1 , H ) : OFFB (H) =A ( 1 , H) : 
OFFE(H)=B(l,H) 

4230  IF  PS (H) =1333  THEN 

MAXB ( H ) =A ( 2 , H ) +C ( 2 , H ) +A ( 1 , H ) +C ( 1 , H ) : OFFB (H) =A ( 2 , H) : 
OFFE (H) =B ( 2 , H) 

4240  IF  PS (H) =1334  THEN 

MAXB ( H ) =A ( 2 , H ) +C ( 2 , H ) +A ( 1 , H ) +C ( 1 , H ) : OFFB (H) =A ( 2 , H) : 
OFFE (H) =B ( 2 , H) 

4250  IF  PS (H) =1335  THEN 

MAXB ( H ) =A ( 2 , H ) +C ( 2 , H ) +A ( 1 , H ) +C ( 1 , H ) : OFFB (H) =A ( 2 , H) : 
OFFE (H) =B ( 2 , H) 

4260  IF  PS (H) =1344  THEN 

MAXB ( H ) =A ( 2 , H ) +C ( 2 , H ) +A ( 1 , H ) +C ( 1 , H ) : OFFB (H) =A ( 2 , H) : 
OFFE (H) =B ( 1 , H) 

4270  IF  PS (H) =1345  THEN 

MAXB (H) =A ( 2 , H) +C ( 2 , H) +A ( 1 , H) +C ( 1 , H) : OFFB (H) =A ( 2 , H) : 
OFFE (H) =B ( 1 , H) 

4280  IF  PS (H) =1355  THEN 

MAXB ( H ) =A ( 2 , H ) +C ( 2 , H ) + A (1,H)+C(1,H) : OFFB (H) =A ( 2 , H) : 
OFFE (H) =B ( 1 , H) 

4290  IF  PS (H) =1444  THEN 

MAXB ( H ) =A ( 1 , H ) +B ( 1 , H ) +C ( 1 , H ) +C ( 2 , H ) : OFFB (H) =B ( 1 , H)  : 
OFFE (H) =A ( 2 , H) 

4300  IF  PS (H) =1445  THEN 

MAXB ( H ) =A ( 1 , H ) +B ( 1 , H ) +C ( 1 , H ) +C ( 2 , H ) : OFFB (H) =B ( 1 , H)  : 
OFFE (H) =A ( 2 , H) 

4310  IF  PS (H) =1455  THEN 

MAXB ( H ) =A ( 1 , H ) +B (1,H)+C(1,H) +C ( 2 , H ) : OFFB (H) =B ( 1 , H) : 
OFFE (H) =A ( 2 , H) 

4320  IF  PS (H) =1555  THEN 

MAXB ( H ) =A ( 2 , H ) +C ( 2 , H ) : OFFB (H) =A ( 2 , H) : OFFE (H) =B ( 2 , H) 
4330  IF  PS (H) >1555  THEN  FA (H) =9 : GOSUB  3080  ELSE  GOTO 
4350 

4340  IF  FA (H) =9  THEN  GOTO  3830 

4350  RETURN 

4360  HW1H=W(1,R,T) 

4370  HX1H=X ( 1 , R, T) 

4380  H Y 1H=Y ( 1 , R , T ) 

4390  HQ1H=Q ( 1 , R, T) 

4400  HZ1H=Z ( 1 , R, T) 

4410  W ( 1 , R , T ) =W ( 2 , R , T ) 

4420  X ( 1 , R , T ) =X ( 2 , R , T ) 

4430  Y ( 1 , R , T ) = Y ( 2 , R , T ) 

4440  Q ( 1 , R , T ) =Q ( 2 , R , T ) 

4450  Z(1,R,T)=Z(2,R,T) 

4460  W(2,R,T)=HW1H 

4470  X (2 , R, T) =HX1H 

4480  Y (2 , R, T) =HY1H 

4490  Q ( 2 , R, T) =HQ1H 

4500  Z (2  , R,  T) =HZ1H 

4510  RETURN 
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